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Chapter 1 
Introduction 

1.1 Rationale for Study 

The abundance of carbon on earth as well as its occurrence in 
celestial bodies has resulted in research over the past century that has 
determined many of the important physical and chemical properties of the 
C 2 molecule. In addition, the relative ease of exciting C 2 has led to 
the discovery of many of its spectroscopic features and molecular con- 
stants. However, the absolute strengths of the band systems of C 2 are 
virtually unknown. This thesis presents the results of a series of 
experimental tests designed to measure the electronic transition moments 
of the Ballik -Ramsay, Phillips, Swan, Deslandres-d'Azambuja, Fox- 
Herzberg, Mul liken, and Freymark systems of C 2 . 

The importance and usefulness of these measurements are stated in 
the following section. A brief description of the problem and the 
experimental technique chosen for its solution are also included. The 
last section of this chapter is a synopsis of the thesis. 

1.1.1 C 2 Band Spectra Occurrences and Uses 

The band systems of C 2 are easily observed in both absorption and 
emission in a wide variety of sources. These sources vary in sophisti- 
cation from an ordinary candle flame to the bov shock wave formed on a 
multimillion dollar planetary entry probe. Several of the more 
important occurrences of C 2 band spectra are discussed below. 

C 2 band systems are excited in virtually all combustion and 
explosion processes and are often used diagnostically in flame and com- 
bustion spectroscopy. Vibrational and rotational temperature 
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determinations from C2 spectra not only indicate the temperature of the 
luminous source but also serve as a diagnostic of local thermodynamic 
equilibrium. Unknown C2 concentrations in a flame or other source can 
be obtained from measurements of the spectral emission of C2 lines or 
band systems if the absolute strength of the band system is known. 

Band spectra of C2 are easily excited by shock heating gases 
containing carbon compounds. This method produces large concentrations 
of C2 in a local thermodynamic equilibrium environment. Consequently 
the shock tube is often used for spectroscopic investigations of C2. 

An additional source of shock excited C2 spectra is the bow shock wave 
on planetary entry probes. Whiting et al. ( 1973 ) have demonstrated 
the feasibility of determining the composition of planetary atmospheres 
by measuring the emission from the heated atmospheric gases in the region 
between the bow shock wave and the probe body. In this technique, a 
measurement of the emission of band systems in selected spectral regions 
can be used accurately to determine the composition of an unknown 
atmosphere. This is only possible if the electronic transition moment 
of the band system being measured is well established. 

Carbon-containing compounds are candidate materials for the heat 
shield of outer planet entry probes. As a result of the high heating 
rates encountered by the probe, the outer surface of the heat shield 
material vaporizes and collects in the boundary layer. These ablation 
products in the boundary layer, of which C2 is a major constituent, 
absorb radiation directed toward the body by the shock layer and can 
substantially reduce the radiative flux reaching the probe surface. The 
band systems of C2 occur in the spectral region from below 0.2 to over 
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2.0 y and thus they can play a major role in absorbing the radiation in 
this wavelength range directed toward the body. The amount of radia- 
tion absorbed is controlled by the relevant absorption coefficients which 
depend on the absolute strength of the respective band systems. Conse- 
quently, in order to evaluate the carbon heat shield effectively, the 
electronic transition moments of the C 2 systems must be accurately 
known. 

Astrophysically, the C 2 molecule is very abundant and of great 
importance (Herzberg, 1950) . The band systems of C 2 are so prominent 
in the atmospheres of the cooler type R and N stars that they are also 
known as carbon stars. Observed band intensities of the C 2 transitions 
can be used to obtain stellar abundances (Fujita, 1970). The intensity 
ratio of the isotopic bands C 12 C* 2 and C 12 C 13 can be used to infer 
abundance ratios of isotopes in stellar atmospheres (Climenhaga, 1960). 

The relative intensities of individual rotational lines in a given 

I 

vibrational sequence and different vibrational sequences of the same / 

band have been used to determine temperatures of stellar atmospheres 
(Aller, 1964). 

C 2 band systems have also been identified in the spectra of the 
photosphere of the sun, sunspots, and comet heads. Therefore, these 
band systems can be used as valuable diagnostic aids by astronomers and 
astrophysicists if the molecular constants are well documented. 

Additional sources of C 2 spectra are discharge tubes containing CO 
and rare gases, electric furnaces, and carbon arc discharges. Phillips 
and Davis (1970) for example have used all three of these sources to 
produce spectra for an extensive analysis of 36 bands of the Swan 
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system. Similar extensive analyses of the other C 2 band systems are 
unfortunately not available. 

Thus it is apparent that there is sufficient interest and need to 
justify a comprehensive experimental program to measure the absolute 
strengths of the band systems of the C 2 molecule. The program to 
accomplish this is described in detail below. 

1.1.2 Definition of Problem 

Simply stated, the goal of this thesis was to measure the absolute 
value of the electronic transition moments of C 2 band systems in the 
spectral region between 0.2 and 1.2 y. Since several of these band 
systems are relatively weak, a source which supplies a high concentra- 
tion (approximately 10 17 particles/cm 3 ) of C 2 particles was required. 

An additional requirement (for ease of data reduction) is that the C 2 
particles exist in a local thermodynamic equilibrium environment whose 
physical properties are well known. A shock tube is one source which 
fulfills these requirements and was selected because of the author's 
familiarity with shock-tube facilities and of the existence of a fine 
spectroscopic shock-tube laboratory at Ames Research Center. Gaseous 
C 2 was produced by shock heating a mixture of argon and acetylene (C 2 H 2 ). 
The radiation from this heated gas was measured using spectrally and 
absolutely calibrated radiometers. The electronic transition moments of 
the band systems were obtained by the synthetic spectrum analysis of the 
absolute intensity data. 

This project is a cooperative effort between Ames Research Center 
and York University. The experimental tests were supported by and done 
at Ames Research Center. TMs work continues the effort of R. W. 
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Nicholls and his colleagues and students at The University of Western 
Ontario and York University to determine basic spectroscopic radiation 
constants of selected molecules. 

1.2 Synopsis of Thesis 

1.2.1 Overview 

Chapter 1 contains the introductory remarks. A synopsis of the 
current knowledge of the C 2 molecule including a brief review of the 
historical background of the band systems is found in Chapter 2. 

The first five sections of Chapter 3 review the parts of shock- 
tube, thermochemistry, and equilibrium radiation theory necessary for 
the understanding of the experimental technique. The final section of 
Chapter 3 describes the synthetic spectrum analysis data reduction 
method. 

Chapter 4 contains a description of the shock-tube facility and its 
operation. The instrumentation and calibration procedures required to 
make absolute intensity measurements of the spectral emission are also 
discussed there. Chapter 4 concludes by presenting a sample of the 
experimental data. 

The measurements of the electronic transition moments of the 
Ballik-Ramsay, Phillips, Swan, Deslandres-d'Azambuja, Fox-Herzberg, 
Mulliken, and Freymark systems are given in Chapter 3. These measure- 
ments are compared with the previous results in Chapter 6. The latter 
part of Chapter 6 summarizes the results and presents final conclusions. 

1.2.2 Sections Reporting New Results 

The measurements of the electronic transition moments are presented 
in section 5.2. An error analysis leading to an estimate of the 
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accuracy of these results is given in section S.3. A critical compari- 
son of the results in section 5.2 with prior work is contained in 
section 6.1. 

For the reader who desires a closer look at the results, Appendix E 
contains tables of the data for each experimental test. These tables 
also include the pertinent thermt chemical quantities for each run. The 
spectroscopic constants and Franck-Condon factors used in the data 
reduction are given in section 2.1 and Appendix C, respectively. 
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Chapter 2 

Synopsis of the Current Knowledge of the C? Molecule 
As mentioned earlier, a vast amount of research and interest has 
been focused on the C 2 molecule. Consequently, many of the details of 
the spectrum including many of the spectroscopic constants have been 
reasonably well determined. This chapter gives a brief historical review 
of the current knowledge of the spectrum of the C 2 molecule. 

2.1 Historical Background of C 2 Band Systems 

C 2 was one of the first molecules to be studied spectroscopically. 
The first band system of this molecule was observed by Wollaston in 1802 
in the spectrum from the base of a candle flame. Wollaston, however, 
was unable to identify the correct emitter. Consequently, the source 
of this band system, now known as the Swan system, remained a mystery 
for many years. Swan (1857), for whom the band was named, observed this 
same band system but concluded that the emitter must be acetylene since 
he could only obtain it by using hydrocarbon flames as sources. Shortly 
after Swan's hypothesis. Van der Willigen (1859) also observed this 
spectrum and properly identified C 2 as the emitter. In the years fol- 
lowing Van der Willigen's work, several other workers agreed with his 
identification of C 2 as the emitter. Among his supporters were 
Schuster (1880), Deslandres (1888), and Kayser (1910). The matter was 
far from resolved however because authorities such as Angstrom and 
ThalSn (1875) and Liveing and Dewar (1880) supported Swan's assertion 
that the emitter was a hydrocarbon. Smithells (1901) suggested that 
neither C 2 nor hydrocarbon were involved but that the proper emitter was 
actually CO. The controversy continued until 1930 when King and Birge, 
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128 years after the first observation of this band system, showed con- 
clusively that the C 2 molecule was the emitter. This demonstration 
resulted from the identification of the C 12 C 13 Swan bands with the 
correct isotope shift clearly resolving the long-term dispute. 

Since the Swan bands were so easily observed in absorption, it was 
assumed that the lower state in the electronic transition was the ground 
state of the molecule. This assumption which was shown to be false 
several years later led to a rather confusing designation of the states 
of C 2 . This problem too was resolved with time as will be discussed 
below. 

While the controversy surrounding the nature of the emitter of the 
Swan bands was still in progress, Deslandres and d'Azambuja (1905) 
observed the second of the band systems of C 2 which now bears their name. 
Mulliken (1930) tentatively identified a third band system of C 2 in the 
ultraviolet. This identification was established in 1939 when both 
Mulliken and Landsverk properly identified the electronic states 
involved in the transition. This system was named for Mulliken. 

Fox and Herzberg (1937) analyzed the fourth band system of the C 2 
molecule. This band system has the same lower state as the Swan bands 
and was found to be strongly overlapped by the nearby Deslandres- 
d'Azambuja and Swan bands. 

Phillips (1948), who perhaps has done more than anyone else to 
increase our knowledge of the C 2 molecule, discovered the fifth band 
system. The Phillips system as it is now known arises from a transi- 
tion between the previously known lower states of the Deslandres- 
d'Azambuja and Mulliken systems. 
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Freymark (1951) observed the sixth C 2 band system. The upper level 
of this transition was a new electronic state which at its discovery was 
the highest lying state of C 2 yet discovered. The lower state is also 
the lower state of the Desl andres-d ' Azambuj a transition. 

Ballik-Ramsay (1959) not only discovered a new band system, the 
seventh, but they showed that the lower state of the Swan transition 
previously assumed to be the ground state was not the ground state of 
C 2 at all. They discovered a new electronic state which lies 610 cm -1 
below the lower state of the Swan transition. The Ballik-Ramsay system 
shares its lower state with the Swan and Fox-Herzberg systems. 

The startling discovery of a "new" ground state by Ballik and 
Ramsay in addition to the new electronic state which is the upper state 
of the Ballik-Ramsay transition led to an inconsistency in the notation 
of the electronic states. Since the notation originally proposed by 
Herzberg (1950) was well accepted, Ballik and Ramsay were reluctant to 
propose a new notation. The situation became more confusing when in 
1963 Ballik and Ramsay discovered a third electronic state lying between 
the upper and lower states of the Swan transition. 

In 1967, Messerle and Krauss discovered the eighth band system. 

The lower state was the same as the Deslandres-d'Azambuja and Freymark 
transitions but the upper state was a new state which arose from 
perturbations on the Oeslandres-d'Azambuja upper state. 

Herzberg, Lagerqvist, and Malmberg (1969) identified three more 
electronic states which lie above all of the other previously identified 
states. The transitions from these new states to lower lying states are 
presently unnamed. With the discovery of these new states, Herzberg 
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et al. proposed a new notation for the electronic states of C 2 . This 
new notation which is now widely in use will be used in this thesis. 
Table 2.1 in the next section which deals with the electronic states and 
their molecular constants contains both the old and new notations for 
reference. 

2.2 Electronic States and Molecular Constants 

A summary of the molecular constants of all of the known electronic 
states of C 2 was published by Ballik and Ramsay in 1963. Since then, a 
new state near the upper state of the Deslandres-d'Azambuja system and 
arising from perturbations on it has been observed by Messerle and 
Krauss (1967). Three more states lying above all of the previously 
identified states were documented by Herzberg, Lagerqvist, and Malmberg 
in 1969. Therefore at the present time there are 14 known electronic 
states of C 2 . The various molecular constants of these states are given 
in table 2.1. A complete listing and bibliography of the data contained 
in table 2.1 are given in the book Spectroscopic Data Relative to 
Diatomic Molecules (1970) edited by Rosen. 

Both the old notation (Herzberg, 1950) and the new notation 
(Herzberg et al., 1969) are given in table 2.1. A perceptive reader 
will notice that the letter B is missing from the singlet notation. 

This letter was reserved by Herzberg et al. for two predicted low 
lying singlet states that have not yet been observed. 

The potential energy curves of the C 2 molecule are shown on 
figure 2.1. All of the electronic states listed in table 2.1 are 
included. The curves represented by the solid lines were constructed 
on the basis of Morse potential functions and were taken from 




Table 2.1 • Molecular Constants for the Known Electronic States of the C 2 Molecule 


State 

T e* 

Ui * 

e 

1 (IX * 

e e 

o) y * 
e 7 e 

B * 
e 

o * 
e 

V 

r * 
e 

Reference 

Old 

notation 

E l Z* 
8 

55034.6 

1671.5 

40.02 

0.248 

1 . 7930 

0.0421 


1.2517 

Freymark (1951) 

• **; 


43240.23 

1829.57 

13.97 


1.8334 

.0204 


1 . 2378 

Herzberg (1950) 

a 't* u 


40796.65 

1106.56 

39.26 

2.805 

1.1922 

.0242 


1.53S0 

Herzberg (1950) 

B 3 *8 

cl *g 

34261.9 

1809.1 

15.81 

-4.02 

1 . 783** 

.0180 


1.2552** 

Herzberg (1950) 

C l% g 

d % 

20022.50 

1788.22 

16.440 

- . 5067 

1.7527 

.01608 

-0.001274 

1 . 2660 

Herzberg (1950) 

A 3 *g 

c K 

13312.1 

1961.6 

13.65 


1.87 



1.23 

Ballik 6 Ramsay (1963b) 

A” 3 r* 

A l w u 

8391.00 

1608. 3S 

12.078 

-.010 

1.61634 

.01686 

-.00005 

1.31843 

Ballik 6 Ramsay (1963b) 

b U u 

b3E 8 

6434.27 

1470.45 

11.19 

.02 

1.49852 

.01634 


1.36928 

Ballik 6 Ramsay (1963a) 

A’ 3 * u 

* 3 * u 

716.24 

1641.35 

11.67 


1.63246 

.01661 

j 

1.31190 

Ballik 6 Ramsay (1963a) 

X' 3 w u 

x *t; 

0 

1854.71 

13.340 

-.172 

1.81984 

.01765 

| -.00023 

1.2425 

Ballik $ Ramsay (1963b) 

! * % 

f l »U 

74550. 

1557. 5 + 



! 1.645 

.019 


1.307 

J Herzberg et al. (1969) 


* \ 

72390. 

1458. 06 + 



1.5238 

.0170 


1.358 

Herzberg et al. (1969) 



70924. 

1360.5 

14.8 


1.448 

.040 ++ 


1.393 

Herzberg et al. (1969) 



*AU units are cm’ 1 except for r e , which are given in \ units. To obtain the T e values for the triplet states, i* 
has £een assumed that the "zero" level of the a % state lies exactly 610.00 cm- 1 above the zero level of the 
X state (Ballik and Ramsay, 1963a). 

**HerzbeTg (1950) quotes different values, which are apparently in error. 

"These are AC (1/2) values. 

* 1.448 - 0. 040fv ♦ 0/2)] ♦ 0.006[v ♦ ( 1/2) J 2 . 
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Figure 2.1.- Potential energy curves for the C 2 molecule. The curves given 
by the solid lines were constructed on the basis of Morse functions. 

The dashed curves were taken from Herzberg et al. (1969) and the dotted 
curve was taken from Messerle and Krauss (1967). 
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Ballik and Ramsay (1963). The dashed curve representing the C' l ir g state 
which results from perturbations of the C l tr g state was adopted from 
Messerle and Krauss (1967). The partial potential curves indicated by 
dashed lines represent the three states observed by Herzberg et al. 
(1969). These curves were drawn on the basis of the data provided by 
them. 

The correct value of the dissociation energy of C 2 has remained 
uncertain for many years. An accurate knowledge of this energy is 
required to predict the thermodynamic properties and C 2 concentrations 
of gas mixtures in which C 2 is a major constituent. It is thus needed 
for the thermochemical calculations which lead to the number densities 
of C 2 in our shock-heated source. Herzberg (1946) suggested a dissoci- 
ation energy value of 3.6 ev. He argued that the high pressure Swan 
bands which arise only from the v ■ 6 level of the upper state of the 
Swan transition were due to an inverse predissociation and that this 
level must be at the same energy as the separated carbon atoms. This 
reasoning is now known to be incorrect (Ballik and Ramsay, 1963). 

Brewer, Gilles, and Jenkins (1948) reported a value of 4.7 ± 0.3 ev 
for the dissociation energy. This value was deduced from measured values 
of the emission of the Swan system at various temperatures using a King 
furnace as a source. Another study by Brewer, Hicks, and Krikorian 
(1962) using this same technique as well as mass spectrometer measure- 
ments yielded an average value of 6.25 ± 0.2 ev. Diowart, Bums, 

DeMaria, and Ingham (1959) performed a mass spectrometric study of the 
gas effusing from a Knudsen cell and listed a value of 6.21 ± 0.07 ev. 
Messerle and Krauss (1967) have made Birge-Spooner extrapolations of 
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the vibrational levels of the d 3 irg and e 3 *g states to limits of 
6.07 and 6.09 ev, respectively. The extrapolation of the E *lg level 
to its limit by Gaydon (1968) indicates a value of 6.4 ± 0.3 ev. 

From these results, it appears that Brewer et al.'s value of 
6.25 ev is reasonable and since Gaydon (1968) also recommends this value 
(6.25 ± 0.20 ev), it will be used for the dissociation energy of C 2 in 
this work. A value of 6.25 ev is shown on figure 2.1 as the dissociation 
energy and therefore represents a slight modification of the curves as 
given by Ballik and Ramsay. 

Theoretical calculations of the electronic structure of C 2 have 
been performed by Mulliken (1939), Clementi and Pitzer (1960), and 
Fougere and Nesbet (1966). The most comprehensive and sophisticated of 
these is the work by Fougere and Nesbet. They have calculated the C 2 
electronic structure over a wide range of internuclear separations using 
three different Slater-type orbital basis sets and two different levels 
of configuration interaction. The calculated values of potential curves 
and spectroscopic constants of 27 electronic states are reported. A 
comparison of these calculations with the experimentally determined 
values given in table 2.1 indicates only a reasonable agreement between 
the two. However, one of the important features of their work is the 
prediction of electronic spectra which contain transitions between pairs 
of states in which only one of the states has been experimentally 
observed. These predictions should be of considerable value in future 
attempts to locate new C 2 spectra. 
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2.3 Electronic Transitions and Spectra 

All of the known electronic transitions of the C 2 molecule are 
listed on table 2.2. The old transition notation when available is 
given in parentheses for reference. Figure 2.2 is a partial energy level 
diagram showing the eight electronic transitions which are of interest 
in this work. The singlet and triplet transitions have been separated 
for clarity. 

Table 2.2 also contains the approximate observed spectral extent 
of each system. Moderate reciprocal dispersion survey spectra showing 
portions of the Phillips, Swan, Deslandres-d'Azambuja, Fox-Herzberg, 

Mul liken, and Freymark band systems are given in section 4.4. 

2.4 Previous Work on Cg Electronic Transition Moments 

This section summarizes the results of the previous work on C 2 
electronic transition moments. Both experimental measurements and 
theoretical predictions are included. A comparison of these results 
with those of the present effort may be found in section 6.1. The 
results are represented in terms of the sum of the squares of the 
electronic transition moments, r | R e / ea 0 1 2 » which is often referred to as 
the electron transition moments squared. 

Other than the present work, the only systematic study of the 
absolute strengths of the C 2 band systems was attenpted by Wentink et 
al. (1967) using the laser excitation technique. They encountered and 
discussed several difficult and conplex problems (see section 6.1) 
which led them to question the validity of their results. They have 
concluded that their results usually give only a lower bound to the 
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Table 2.2 Band System Designations of All Known Electronic 
Transitions of the C 2 Molecule. 


Band system 

Electronic transition. 
Upper state - Lower state 

Approximate 
spectral # 
region, A 

Swan 

d 3 v g - a 3 ir u 
(A 3 iTg - X* 3 ir u ) 

4200 - 7000 

Deslandres -d ' Azambu j a 

C ^i?g — A 
(c - b 

3300 - 4100 

Fox-Herzberg 

e \ - a 3 * u 

(B 3 ir g - X' 3 ir u ) 

2300 - 3300 

i 

Mull i ken 

! 

( d 1 £ J . X >£*) 

2200 - 2500 

Freymark 

Eh g' A 
(e - b *ir u ) 

2000 - 2400 

Phillips 

A l w u - X h* 
(b\-X 'tp 

7000 - 12100 

Ballik-Ramsay 

b 3 Eg - a 3 ir u 
(A* - x* 3 * u ) 

11600 - 25000 

Ultraviolet bands observed 
by Herzberg et al. 

(1967) 

f 3 Eg - a 3 ir u 
g 3 A g - a 3 ir u 
P - X lE g 

1300 - 1450 

Messerle-Krauss 

C - * **U 

3100 - 3400 


























- 17 - 


■i 



Figure 2.2.- Energy levels diagram for C 2 showing the inportant electronic 
transitions. The singlet and triplet transitions have been separated 
and some of the electronic states omitted for clarity. 
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electronic transition moments squared value. Consequently on the 
following pages, their results will be given as lower bounds. 

Since the Swan system is the strongest radiating system of the C 2 
molecule, it has been studied extensively. Table 2.3 is a comparison of 

the experimental values of the E|R e /ea 0 | 2 value of the Swan system 

found in the literature. Only measurements of the Mr * 0 sequence 
(5165 X) are included. 

The shock-tube results of Arnold (1968) and Fairbairn (1966) 
obtained behind the incident shock agree well with each other and rea- 
sonably well with the results of Harrington et al. (1966), Sviridov 
et al. (1965), and Sviridov et al. (1966) using a shock tube in the 
reflected shock mode. With the exception of the measurements by Hicks 
(1957) in a King furnace, the data obtained in furnaces appear too small. 
The lifetime measurements of Jeunehomme and Schwenker (1965) and Fink 
and Welge (1967) differ by a factor of 4 with the latter measurement in 

fair agreement with the shock-tube data. The results of Wentink et al. 

(1967) which are also lifetime measurements represent a lower bound. 
Consequently, the measurements of the absolute strength of the Swan 
system differ by a factor of about 6. However, it is generally accepted 
(Arnold, 1968) that the true E|R e /ea 0 | 2 value is probably bounded by 
the shock-tube measurements listed in table 2.3. 

A recent measurement of the variation of the Swan system with inter- 
nuclear separation was reported by Danylewych (1971) and Danylewych and 
Nicholls (1974) . This work contains an excellent summary and comparison 
of all of the past measurements of the variations of spectral intensity 
of the Swan system. 
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Table 2.3 Measurements of the Values of the Electronic Transition 
Moments Squared of the A v ■ 0 Sequence (S165 A) of the C 2 
Swan System. 


Experimenter 

Method 

l|R e /ea 0 | 2 , atomic units 

Arnold (1968) 

Shock tube 

3.56 ± 0.50 

Fairbaim (1966) 

Shock tube 

3.36 ± 1.20 

Harrington et al. (1966) 

Shock tube 

2.95 ± 0.92 

Sviridov et al. (196S) 

Shock tube 

2.24 ± 0.81 

Sviridov et al. (1966) 

Shock tube 

2.85 ± 0.20 

Hicks (1957) 

King furnace 

2.64 

Hagen (1963) 

King furnace 

0.71 ± 0.40 

Jeunehosroe 6 Schwenker (1965) 

Laser excitation 

0.60 i 0.02 

Fink 8 Welge (1967) 

Phase shift 

2.24 ± 0.70 

Wentink et al. (1967) 

Laser excitation 

1.35 
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The results of various theoretical efforts to predict the absolute 
strength of the Swan system are listed in table 2.4. These results show a 
discordance similar to that of the experimental measurements. The early 
predictions of Shull (1951) and Stephenson (1951) show surprisingly good 
agreement with the higher values reported from shock- tube tests. 

Hie experimental and theoretical results of the absolute strengths 
of the other C 2 band systems of interest are summarized in table 2.5. 

It is readily apparent that very little experimental or theoretical work 
has focused on these weaker band systems. Of these six systems, the 
only one whose strength has been measured absolutely is the Mul liken 
system (Smith, 1969). The measurement of the Phillips system by Hicks 
(1957) was made relative to the (1,0) band of the Swan system and has an 
uncertainty of a factor of 3. The measurement of the Fox-Herzberg sys- 
tem by Sviridov et al. (1969) is an upper limit. As discussed earlier, 
the measurements of Wentink et al. (1967) represent lower bounds. 

Clement! (1960) using the dipole-moment operator method calculated 
the absolute strength of five of these band systems. His results are 
presented in the form of oscillator strengths or f numbers. It should 
be noted that Clement! apparently omitted the statistical weight of the 
upper state from the definition of the f number. Consequently, the 
values given in table 2.5 have been modified to account for this factor. 
With the exception of the work on the Deslandres-d’Azambuja system by 
Shull (1951) and Stephenson (1951), no other theoretical calculations of 
the strength of these systems are known to the author. 

Most of the results summarized in tables 2.3 to 2.5 are reported in 


terms of oscillator strengths or f numbers. The astute student of 
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Table 2.4 Theoretical Predictions of the Absolute Strength of 

the C 2 Swan System. 


Reference 

Z|R e /ea 0 | 2 . A.U. 

Remarks 

Clementi (1960) 

1.14 

E|R e /ea 0 | 2 value modi- 
fied to account for the 
statistical weight of 
the upper level omitted 
from defining equation 

Shull (19S1) 

4.24 

Same as above 

Stephenson (1951) 

4.06 


Lyddane, Rogers 5 Roach 
(1941) 

3.39 


Spindler (1965) 

5.04 

Based on Jeunehomme*s 
lifetime measurement on 
the d state 

Coulson 5 Lester (1955) 

33.92 



Table 2.5 Experimental and Theoretical Results of the Absolute 
Strength of the Band Systems (Excluding the Swan System) of 
the C 2 Molecule. 


Band system 

Results in literature *E|R e /ea 0 | 2 , A.U. 

Experimental 

Theoretical 

Ballik-Ramsay 


1.21 - Clementi (I960)* 

Phillips 

0.68 - Hicks (1957) 

0.13 - Clementi (I960)* 

Deslandres- 

d’Azambuja 

> 0.18 - Wentink et al. (1967) 

1.12 - Clementi (I960)* 
3.97 - Shull (1951)* 

1.35 - Stephenson (1951) 

Fox-Herzberg 

< 0.55 - Sviridov et al. (1965) 
> 0.45 - Wentink et al. (1967) 

57.80 - Clementi (I960)* 

Mulliken 

0.42 - Smith (1969) 

0.78 - Clementi (I960)* 

Freymark 

— ni .J 




*These values are modified to account for the statistical weight 
of the upper level omitted from the defining equation. 
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spectroscopy realizes that this nomenclature is ambiguous and that 
caution must be exercised when using results described by this terminol- 
ogy. Consequently the author has investigated the original references 
listed on these tables to determine the manner in which the reported 
f number was defined. Therefore the results originally reported as 
f numbers have been converted to E|R e /ea 0 | 2 values by using a consistent 
set of defining equations. These relations were taken from Schadee 
(1967) and are given in Appendix A. 

2.5 Goals of the Present Work 

It is rather obvious from examining table 2.5 that there is a real 
need for a systematic determination of the absolute strength of the band 
systems of the C 2 molecule. Thus the present effort was motivated by 
the gross uncertainties that exist in the absolute strengths of the 
band systems of C 2 . The previous section clearly demonstrated this 
uncertainty. 

The goal of this work is to obtain measurements of the absolute 
strength of the band systems of the C 2 molecule. These results which 
are reported in terms of electronic transition moments squared values 
are deduced from absolute measurements of the radiation by the synthetic 
spectrum analysis technique. The interest and usefulness of these 
results was previously set forth in section 1.1. 


j> 
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Chapter 3 

Theoretical Bases of Experiments 

This chapter is a review of the theory which is necessary for the 

understanding and interpretation of the experimental measurements. The 

most important parameter in this work is the integrated intensity, E^ u , 

from spontaneous emission in a rotational line emanating from a light 

source in local thermodynamic equilibrium. The general equation for 

this intensity is (Nicholls, 1969): 

16 x 10 -7 

E fc U = 3(2 - 6 77T"C2S* '♦ 1)(2J* + D^iv'A’J' * Z l R e (r v'v"^ 2 * Vv" 

J'A' 

• S (3.1) 

j"A" 

The principal goal of this work is to determine the quantity 

z|R e (f v » v '0 I 2 which is the square of the electronic transition moment of 

the band system. This is accomplished as follows: Firstly, the absolute 

value of the intensities E£ U of several band lines are measured using 

calibrated radiometers; secondly, the quantity, N uv » a »j», which is the 

number of excited emitters involved in the radiative transition is 

obtained from the thermodynamic properties of the source (in this work, 

the shock-heated gas behind an incident shock wave); thirdly, the neces- 

J'A' 

sary Franck-Condon factors, q v , and HOnl-London, S , factors are 

v J"A" 

calculated from the theory of molecular spectroscopy. The remaining 
parameters in equation (3.1) can be obtained fr v >m the known molecular 
constants of C 2 . Thus the quantity £ | R (r, ,,,„)| 2 can then be deter- 
mined directly from equation (3.1) after summing it over all of the lines 
which contribute to the measured radiation. 
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The structure of this multicomponent chapter is as follows: 

Equation (3.1) is derived from equilibrium radiation theory in section 

3.1. The calculation of the Franck-Condon and HBnl-London factors is 

* 

discussed in section 3.2; and section 3.3 reviews the theory of radiative 
transfer. The application of the general equation (3.1) to the specific 
transitions and band systems of C 2 is reviewed in section 3.4. Section 
3.5 reviews the parts of shock tube and thermochemistry theory necessary 
to provide the equilibrium concentrations n uv «^ij» as we ^ 1 as t ^ e 
equilibrium temperature of the source. The final section (3.6) of this 
chapter reviews the synthetic spectrum data analysis technique and 
relates the theory of the preceding five sections to the analysis of the 
experimental data. 

3.1 Equilibrium Radiation Theory 

The principles of equilibrium radiation theory that pertain to 
diatomic molecules are reviewed here. Full discussions of the quantum 
mechanical model of radiative transitions in a diatomic molecule have 
been made by Pauling and Wilson (1935), Herzberg (19S0), and Nicholls 
(1969) . 

3.1.1 The Molecular Wave Equation 

The theoretical study of radiative transitions in diatomic molecules 
involves the formulation and solution of the molecular wave equation. 
Schrbdinger's time independent equation for a system of particles repre- 
senting a diatomic molecule is 

HY » EY (3.2) 

where H is the Hamiltonian operator, Y is the total wave function, 
and E is the eigenvalue which represents the total energy of the 
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system. Hie solution of equation (3.2) yields a series of Y's and E's 

from which the quantities E|R (? , ,,) | 2 , q v , v „, v » and S J A in 

“ v v v v uXr J m A h 

equation (3.1) can be calculated. The Hamiltonian is the total energy 
of the system and is usually approximated by the sum of the kinetic and 
potential energies of the nuclei and electrons, that is, 

H - (K.E.J n . (K.E.) C . V„ . V m . V M 


■-$( 


.ybh^L. yb' 2 

T jy r 


(3.3) 




n n,e 

The potential energy terms represent the electron- electron ( v ee )» 

nucleus-nucleus (V ), and electron-nucleus (V ) Coulombic interactions, 
nn ne 

In equation (3.3), the mass of the nuclei and electrons are represented 

by M and m , respectively, while the quantities r.. f , r.. f , and r. . 
n c lijjij 

are the distances between the pairs of interacting particles. The 

Schrbdinger wave equation for this multiparticle system now becomes 
X V _ 2 — « V 2 


y*i.y±.,.± 


-ft 2 


(V + V + V )V » EV 


ee 


nn 


ne' 


(3.4) 


n e 

The Born-Oppenheimer approximation which asserts that the true 
molecular wave function may be represented by the product of two differ- 
ent wave functions - one for the electronic motion and the other for the 
nuclear motion, that is. 


’total < r - O ' < 3 - 5 > 

is used to obtain solutions to equation (3.4). The symbol K represents 

the electron coordinates and r is the internuclear separation. Also 
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inherent in the Born-Oppenheimer approximation is the assumption that 
the total energy can be separated into two distinct parts, that is, 


E , a E ♦ E . 
total e n 


(3.6) 


Substituting these approximations into equation (3.4), expanding the 
Laplace operators, dividing by the total wave function, and neglecting 
terms of the order of m/M gives an equation whose terms can be separated 
into two wave equations. One equation describes the electron motion, 

£V*. (r ' 5) * ff \' E e (r) - v ee (r - 0 ' v „„< r - «lV r < » ' »• 


and the other describes the nucleav motion, 

. ill [En(r , . V nll (r)J, n Cr) - 0. 


(3.7) 


(3.8) 


n 


The quantities in brackets in equations (3.7) and (3.8) are usually 
rewritten in terms of the molecular potential U(r) which is the sum of 
the electronic eigenvalue E e (r) and the nuclear potential V nn (r). The 


separated wave equations are now given by 

E V e 2 V» e (r, O ♦ ^~® e [U(r) - V(r)]i|> e (r, £) * 0, 

n ^ 


and 


n 


r_( r ) e-2 

+ — [E - U(r)]* (r) = 0, 
h 2 " 


Mn 


(3.9) 


(3.10) 


for the electron and nuclear motions, respectively, where V(r) is the 
sum of the potential energy terms. 

The solution of equation (3.9) is a very difficult problem. Equa< 


tion (3.10) is somewhat easier to solve although it involves a knowledge 
of the molecular potential function U(r) which is obtained theoretically 
from the solution of equation (3.9). However, at the present time. 
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theoretical calculations of U(r) exist for only the simplest molecules. 
Therefore, the normal procedure is to adopt an approximate form for U(r) 
which contains parameters that are determined by a comparison of observed 
and calculated energy levels and use this analytical form to assist in 
the solution of equation (3.10). 

With the assumed form of the electronic energy, equation (3.10) is 
solved using various methods and models of varying sophistication to 
obtain the vibrational and rotational energy levels. In this study, we 
shall use the energy eigenvalues of the vibrating-rotator model. The 
vibrational and rotational energy values given by Herzberg (1950) are. 


respectively. 


hc [“e( v * l) ■ Ve( v * t) ! * Ve( v * * Ve( v * l)\ 


hcG(v) , 


Er * hc[B v (J + 1) - D v J 2 (J ♦ l) 2 + (A - B V )A 2 ] 
* hcF(J) , 


(3.11) 


(3.12) 


where 


B = B 
v e 


- ° e ( v * i) 


B e - — 

8ff 2 cl, 


D = D + B„ 
v e e 


M) 


In the above equation, w t w x , <*> y , and <*• z are vibrational constants, 

0 G 6 6 0 6 6 

B, D, and ot are rotational constants, and the subscript e refers to 
0 

the equilibrium position. 
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The total energy of the molecule can therefore be written as 

E total * hc[T e + (3.13) 

where T e is the energy of the minimum of the electronic potential well 

of the electronic state of interest. This equation is used in the fol- 
lowing section to calculate the central wavelength of the rotational 
line whose total emissive intensity is given by equation (3.1). 

3.1.2 Electronic Transition Wavelengths 

When a transition from a upper electronic state to a lower electronic 
state whose energies are given by E' and E" takes place, a photon of 
energy hv is emitted (Bohr’s frequency condition), that is, 

E' - E” = hv = hc/X (3.14) 

In this equation, v and X are the frequency and wavelength of the 
emitted radiation. From this relationship and the molecular energies of 
the previous section, the central wavenumber of a spectral line resulting 
from a transition between two electronic states is 

v(cnT l ) - (T ' - T ") + [G’(v') - G"(v”)] ♦ [F’(J') - F’’(J'')]. 

(3.15) 

This quantity is one of the input parameters required for the 

successful evaluation of equation (3.1). 

Equation (3.15) is often written as 

v(cm _1 ) » v (v’, v") + [F'(J') - F"(J")] - (3.16) 

° X(A) 

where v q (v’, v") is called the band origin of the (v', v") transition. 

Calculated band origins (in vacuo) for the (v', v”) transitions of the 

C 2 molecular systems of interest in this study are listed in appendix B. 



-29- 


In order to calculate the central wavenumber of a rotational 
line (equation (3.15) or (3.16)), an accurate expression for F(J) must 
be used. The rotational term energies, F(J), depend on the manner in 
which the various angular momenta couple together and are discussed in 
the next section. 

3.1.3 Hund's Coupling Cases 

The interactions of the various angular momenta are classified on 
the basis of the strength of the coupling to the intemuclear axis and 
are referred to as Hund's coupling cases (a), (b), (c), (d), and (e). 

Since the C 2 transitions of interest in this work are represented by 
either Hund's case (a) or (b), only these two cases will be discussed. 
Herzberg (1950) gives a more detailed discussion of all of Hund's cases. 

The vector notation of H».rzberg (1950) is used in the following discussion. 

Hund's Case (a) 

In this case, it is assumed that the coupling of the nuclear rota- 
tion with the electronic motion is weak compared to the coupling of the 
electronic motion to the intemuclear axis. Figure 3.1 is a vector 
diagram for Hund's case (a). As shown in this figure, the rapid pre- 
cession of the orbital angular momentum t and the spin angular 
momentum $ about the intemuclear axis result in the axial components 
t and ft. The total electronic angular momentum along the intemuclear 

4 + 

axis is therefore SI * A + E whose values vary from A - E to A ♦ l. 

The total angular momentum ft results from the coupling of ft with 
the nuclear rotation vector ft. The latter two vectors nutate about ft 
as indicated in the figure. The total angular momentum is constant in 
magnitude and direction and always exceeds its component ft. This 
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Figure 3.1.- Hund's case (a). 

accounts for the fact that energy levels with J < 0 are not observed. 
The rotational energy eigenvalues for this case now become 

E r - hc[B y J(J ♦ 1) - D v J 2 (J ♦ l) 2 - (A - ByJO 2 } (3.17) 

where 

J • fl, ft ♦ 1, fl + 2, . . . 

Hund's Case (b) 

Hund's case (b) is characterized by a very weak or zero coupling of 
the spin angular momentum to the internuclear axis. Therefore only t 
couples to the internuclear axis to produce its component X along the 
internuclear axis. This vector, t, and $ couple to form an intermediate 
resultant it which is called the total angular momentum apart from spin. 
The spin angular momentum, §, now couples to £ to form the total 
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angular momentum vector, J. A vector diagram of this coupling case 
illustrating the vector motions is shown in figure 3.2. As a result of 
the coupling between £ and a splitting of the energy levels with 
different J but equal K occurs. This phenomenon is observed in the 
rotational levels of the Swan, Fox-Herzberg, and Ballik-Ramsay systems. 
Consequently the rotational energy eigenvalues take different forms for 
the various transitions and will be discussed later when the theory is 
applied to the specific C 2 transitions. 

It should also be noted that Hund's coupling cases are limiting 
cases and quite often as the rotation increases, a transition from one 
coupling case to another takes place. This is because the increased 
rotation rate tends to uncouple the angular momentum vectors from the 
intemuclear axis. 



Figure 3.2.- Hund's case (b). 
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3.1.4 Electronic Transition Intensities 

The wavelength of an emitted spectral line resulting from a 
transition between electronic states was obtained in the previous sec- 
tions. The probability that this transition will take place, which 
leads to the derivation of equation (3.1), is discussed below. 

The Einstein transition probability of spontaneous emission, A^, 
which is the total probability per second per particle that a spontane- 
ous transition from an upper state u to a lower state A will occur 
accompanied by emitted radiation of frequency * [(E u - E^)/h] 
is given by (Herzberg, 1950), 

A ui " — IH ” 1 l^l*. (3.18) 

If one considers only electric dipole transitions, then the matrix 
element R which is known as the transition moment is defined as 

S 0 * * / ¥ u *frr £ dx (3.19) 

ii £ 

where ¥ and ¥ are the total eigenfunctions of the states involved in 
the transition and M is the electric dipole moment given by 

- £ er r (3.20) 

i 

By writing the total eigenfunctions as the product of electronic, 
vibrational, and rotational components and transforming from electron 
coordinates to molecular fixed coordinates, Nicholls and Stewart (1962) 
showed that the transition moment for a molecular line could be written 
as 

R iv"J ,, A"m" | * er i|*t)*v" dr 

* /^‘A'm 1 ^*J ,, A H m" sin ® d ® d ®* 
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Squaring this equation and sunning over the appropriate states m' and 


m" one has 


I „uv'J , A , l l 2 1 1 » 1 2 

R tv"J"A"| * |J Wv" df | * s 


^ * s independent of r , it can be removed from the integral and 
equation (3.22) becomes: 


jUv'J'A' 2 

tv"J"A" 


l R e 1 2 ’ |/W *| 2 * S 


(3.23) 


However, R e is usually dependent on r. Since this dependence can 
often be represented by a power series in r, Fraser (1954) defined a 
characteristic internuclear separation, r y , y „ known as the r centroid 


/» v t» v ..r dr 

/ K'K" dr 


(3.24) 


Using this approximation, equation (3.22) becomes 

I uv 1 J 1 A 1 ^ | ^ 12 

R 1v"J"A m * l R ef r V 'v"^ 2 * /♦y'V dr * S . 


l R e^v'v M ^ 2 * %r«v M ® * 

. J"A" 


(3.25) 


where, l R e (* v . v ..) 1 2 is the electronic transition moment squared, r is 

the r centroid, q^ f „ is the Franck -Condon factor, and S’?,'*' is 
the HBnl-London factor. Each of these quantities will be discussed in 
detail in subsequent sections. The product of | R e (r v , v „) | 2 • q^,, is 
often referred to as the band strength. Equation (3.18) for A can 

VI Mt 


be written as 


64v**v* t f 

A ut * Hnff" S l R e (f vV ) l 2 * Vv" * S j"A" * 


(3.26) 


where |R e (^ v# ylf ) | 2 - l\%(i ytyt ,)\ 2 /du. 
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The total energy emitted during a spontaneous emission transition 


from a upper state u to a lower state A is given by 

« u 


10-7 . 

- 3 — hv ,N A ,, watts/cm 3 »sr 
4 ir u£ u u£ 


(3.27) 


where is the number density of emitters in the upper state. Equa- 


tions (3.26) and (3.27) allow the equation for the emitted energy in a 
single line to be written as 
.. u 


16 x 10’W.c 
u£ 


^uv'A'J' E l R e ^v'v"^ l 2q v'v" S J"A" 


3(2 - « 0 >a O(2S' + 1)(2J' ♦ 1) 

(3.28) 

which is identical to equation (3.1). In this equation, the total 
statistical weight is the product of the rotational statistical weight 


(2J* + 1) and the electronic statistical weight (2 - ^,)(2S' ♦ 1 ). 


The Kronecker delta symbol, 6 . suggested by Schadee (1964) is adopted 

and is zero when i t j and unity when i = j. The vibrational statistical 
weight is 1 . 

Since the bandpass of the experimental measurements contained 
several rotational lines, equation (3.28) must be summed to include all 
lines which contribute in the Spectral region of the bandpass. Similarly 
by summing over all states involved in an electronic transition, the 
total emitted energy in that transition due to spontaneous emission is 
given by 


E u , Z 1 * * lo ~ 7 V,c0 4 ) y T — q -y-. v " 

x |_p'p" A ' A " ( " o.A ,)(2S 


_ J'A* ~ 

>A " fN 

2j‘ ♦ 1 1 uv' J'A* • 


(3.29) 
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The 2!) accounts for the parity of the states. It should be noted 
p ' p " J'A» 

that the Httnl- London factors Sj„^„ in equation (3.29) are normalized 
in the manner suggested by Schadee (1967): 

2 E S j»A" * (2S ' + 1)(2J ' + lh (3,30) 

A 'A" J" 


3.2 Franck-Condon and HBnl- London Factors 

In equation (3.25), it was asserted that the band strength could be 
separated into two factors through the r centroid approximation of 
Fraser (1954). That is 


1/ 




- IMW'I 2 ’ (3 - 3I) 


where r was defined as 


/ W r dr 


r, v" P 

/ VV' dr 

and the ip * s are the vibrational wave functions and r is the inter- 
nuclear separation.. Despite some criticisms of applicability to individ- 
ual cases (James. 1966 and Klemsdal, 1973) this approximation has been 
shown to be adequate for bands of most systems (McCallum, Jarmain, and 
Nicholls, 1970 and Nicholls, 1974). 

The Franck-Condon factor 


( 3 . 32 ) 


flylyH 


1/ 


^y I ’fyt 


dr 


(3.33) 


is the square of the overlap integral of the vibrational wavefunctions 
and its magnitude is controlled by the relative position of the poten- 
tials and the relative phase of the wavefunctions. It exerts a very 
strong influence on the relative band intensities within a system. In 
the present work, the Franck-Condon factors of McCallum, Jarmain, and 
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Nicholls (1970) are used. These Franck-Condon factors are based on 
realistic Klein-Dunham potentials derived from measured molecular con- 
stants. A complete list of the Franck-Condon factors used in this work 
is given in appendix C. 

The HUnl -London factors were defined as 


S J"A" = / ' A ' I D I < * , J"A ,,dr 


(3.34) 


where ip j,^, and ar © the electronic-rotational wave functions for 

the upper and lower levels, and D is the dyadic of coordinate trans- 
formation. These factors are well known for most common types of transi- 
tions. The H8nl -London factors used in this work were taken from 
Schadee (1964) whose calculations have been verified by Kovacs (1969) 
and Whiting (1973). 

The electronic transition moment z|R (r , „)| 2 is often 
nondimensionalized by dividing by the product of the electron charge, e, 
and the Bohr radius, a_. The band strength is then written as 


Vv" = E|R e ( Vv” )/ea ol 2 * Vv" = z l R e /ea o * 2 * Vv"* A ’ U ' 


(3.35) 


The total line strength can be defined as 


C’j" = E l R e /ea ol 2 * Vv" * S J A * A - U " 


(3.36) 


where roughly speaking the electronic transition moment, e|r /ea | 2 , 
determines the absolute strength of the whole band system, the Franck- 
Condon factors, control the relative strength of the bands within 

j»A» 

a system, and the Honl- London factors, S JIIAII , 
intensities of the lines within a band. 


determine the relative 
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3.3 Radiative Transport Theory 

The data reduction technique used in this study requires the 
theoretical computation of spectral intensities. In order to compute 
these intensities accurately, it is necessary to consider the interaction 
of the radiation (given by equation (3.1) or (3.28)) with the heated gas 
particles. This interaction is well described by the equation of radia- 
tive transfer and is discussed in detail by Vincenti and Kruger (1965). 

The one-dimensional radiative transfer equation which applies well 
to the case of shock -heated gas is (Vincenti and Kruger, 1965) : 

ar * •»[> - M- iw)]v T) * “A «•>(- nr) - Vx ■ (J - S7) 

where 1^ is the specific intensity along a given path x, B x is the 
Planck black-body function, and <* x is the absorption coefficient. This 
equation is based on the assumption that the gas is in local thermody- 
namic equilibrium and that the various atomic and molecular states are 
populated predominately by particle collisions rather than through radia- 
tive emission or absorption. The three terms on the right-hand side of 
equation (3.37) are due to spontaneous emission, induced emission, and 
absorption, respectively. Integrating equation (3.37) over x and 
substituting c^’ * a x [l - exp(-hc/XkT)] , which is the absorption 
coefficient corrected for stimulated emission, one gets 

I x = B X (T)[1 - exp(-ot x H)] + I x (A s )exp(-a x ')0, (3.38) 

where & is the optical path length, l s represents the gas boundary, 
and Ix(* s ) * s specific intensity incident on the gas boundary. 

If one also assumes, as is true in the present tests, that the 
spontaneous emission processes dominate the absorption and induced 
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emission processes, equation (3.38) becomes 

I x - B X (T)[1 - exp(-o x 'A)]. (3.39) 

In the present tests the gas is also optically thin, that is, a^' « 1, 
and the exponential in equation (3.39) can be expanded to give 

I A = a x ’AB x (T). (3.40) 

Using Kirchoff's law, a x ' = E X /B X (T), this becomes 

I x = E x *. (3.41) 

The spontaneous emission intensities (equation (3.1) or (3.28)) 
described in section 3.1 can now be used to compute the specific inten- 
sity directly. With this information, a synthetic spectrum (section 3.6) 
can be calculated and consequently the data reduction procedure can be 
completed. 

3.4 Application of Theory to C 2 Transitions 

The theory of molecular spectroscopy reviewed above is here applied 
to the specific transitions of C 2 . However, before this application is 
made, the selection rules which limit the transitions between energy 
levels must be recalled. 

3.4.1 Electric Dipole Transition Selection Rules 

Selection rules govern both the quantum numbers and the symmetry 
properties of the states involved in a transition. In the following 
review, the selection rules are valid for both Hund's case (a) and 
case (b), except where noted. 

The symmetry selection rules depend on the symmetry properties of 
the eigenfunctions of the state under consideration. A rotational level 
is designated positive or negative depending on whether or not the total 
eigenfunction changes sign when it undergoes a reflection at the origin. 
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This property is known as parity. A rotational level of a homonuclear 
molecule is also designated symmetric (s) or antisymmetric (a) if when 
the nuclei are interchanged the total eigenfunction remains unchanged or 
only changes its sign. The symmetry selection rules of the rotational 
levels for a dipole transition in a homonuclear molecule are: 

(1) Positive levels combine only with negative levels and vice 
versa, that is, 

+ «-*-, + +, -«/*•-. 

(2) Symmetric states combine only with symmetric states and anti- 
symmetric states combine only with antisymmetric states, that is, 

sym — *■ sym, sym antisym, antisym — *• antisym. 

Since the C 2 molecule has zero nuclear spin, the latter selection rule is 
true for all transitions. 

In addition to the symmetry properties of the total eigenfunction, 
the symmetry properties of the electronic eigenfunction are also important 
in determining the selection rules of electronic transitions. These 
properties depend on the symmetry properties of the field in which the 
electrons move. A nondegenerate state (E state) is called positive (E + ) 
or negative (£~) if its electronic eigenfunction remains unchanged or 
changes sign when reflected at any plane passing through both nuclei. 

For degenerate states, linear combinations of the eigenfunctions can be 
obtained in order that the electronic eigenfunctions sign either changes 
or remains the same. Consequently this nomenclature (+ or -) is usually 
omitted for degenerate states. If an electronic eigenfunction remains 
unchanged when reflected at the center of symmetry (midpoint of the 
internuclear axis for a homonuclear molecule), then the state to which 
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it belongs is called an even or "gerade" state and is indicated by the 
subscript g. If the eigenfunction changes signs when reflected at the 
center of symmetry, its state is called an odd or "ungerade" state and 
is indicated by the subscript u. The selection rules resulting from 
the electronic eigenfunction properties are: 

(1) Positive states combine only with positive states and negative 
states only with negative states, that is, 

I + *-*■ Z* , Z~ E~, Z + +b- Z ~ . 

(2) Even electronic states can only combine with odd states or 
vice versa, that is, 

g «-*■ u, g •*/*• g, u u. 

The selection rule for the total angular momentum number, J, is 

AJ = 0, ±1 

where J = 0-*-J=0is prohibited. The selection rule for the quantum 
number A which represents the component of the electronic orbital 
angular momentum along the internuclear axis is 

AA = 0, ±1. 

An additional selection rule imposed on the resultant spin, S, is 

AS = 0. 

This rule implies that only states of the same multiplicity can combine 
with one another. 

As mentioned above, the previous selection rules apply generally 
to both Hund's case (a) and case (b). The following selection rules 
apply only if both of the electronic states involved in the transition 
belong to the same coupling case. For Hund's case (a), the following 
selection rules are valid: 
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(1) AE « 0 

(2) Aft * 0, ±1 

where £ and ft are the quantum numbers which represent the component of 
spin along the internuclear axis and the component of the total electronic 
angular momentum about ''.he internuclear axis, respectively. If the 
quantum number ft is zero for both electronic states, then an additional 
restriction is placed on the quantum number J - only AJ - ±1 transitions 
can occur. 

In Hund's case (b), the following selection rule is placed on the 
quantum number K, the total angular momentum apart from spin, 

AK = 0, ±1 

where AK - 0 is forbidden for I-E transitions. 

All of these selection rules are summarized in table 3.1. 

The structure of each of the C 2 band systems studied experimentally 
will now be discussed. 

3.4.2 The C 2 Swan System Transition, d 3 ir„ - a 3 ir u 

The C 2 Swan system (Swan, 18S7) results from transitions between 
the d 3 iTg and a 3 ir u electronic states (see figure 2.1). Term energies 
corresponding to these states are given in table 2.1. The vibrational 
term energies are calculated using equation (3.11) and the vibrational 
constants which are also listed in table 2.1. A tabulation of the band 
origins, A 0 , of this system in appendix B shows that the approximate 

o 

spectral range is between 3600 and 8600 A. 

Franck-Condon factors (appendix C) which control the relative 
intensity distribution of the bands within a system can be used to 
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Table 3.1 Electric Dipole Transition Selection Rules. 
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determine the most intense bands. Figure 3.3 is a Deslandres diagram 
of Franck-Condon factors normalized to the maximum value. Only those 
normalized values exceeding 0.1 are included. As observed on the figure, 
the most intense bands form a parabolic curve whose axis is the principal 
diagonal of the diagram. This curve is known as the Condon parabola or 
Condon locus. The width of this parabola increases with increasing 
difference (Ar e ) of the equilibrium internuclear distances (minima of 
the potential curves) of the upper and lower states. Since Ar e is 
small for this system, the narrow width of the Condon parabola for the 
Swan system is expected. Figure 3.3 can be used a priori to determine 
rapidly the prominent features of the excited Swan system. 

The Swan system is an example of a transition in which both states 
are well represented by Hund's case (a) at low values of J and Hund's 
case (b) at large values of J. Figure 3.4 is an energy level diagram 
of the rotational structure of the Swan system. Both electronic states 
are assumed to comply with Hund's case (b) conditions. The rotational 
energy levels, where Fj, F 2 , and F 3 refer to the levels with J ■ K + 1, 
K, K - 1 , respectively, are given by 

Fi = B v K(K ♦ 1) ♦ D v K 2 (K + l ) 2 ♦ r J T ♦ k(y - 2 IV 3 ) + M K ) > 

F 2 = B y K(K ♦ 1 ) ♦ D v K 2 (K + 1)2 - .-y - y + 2 e - $ ± (K) , >(3.42) 

F 3 = B v K(K ♦ 1 ) + DK2(K + 1)2-|- (K + l)(y ♦ ♦ * ± 00 , 

where 4 > + (K) is the A doubling function and is different for the two 

(±) components. In the above equations, A is the coupling constant and 
y and e are splitting constants. There are three possible values of J 
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Figure 3.3.- Deslandres diagram of normalized Franck-Condon factors for 
the Swan system. Only normalized values exceeding 0.1 are included. 
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Figure 3.4.- Energy level diagram for the rotational structure of the 
Swan (d 3 ir- - a 3 ir u ) transition K" * 12. Both electronic states 
are inverted and assumed to belong to Hund's case (b) . The anti- 
symmetric levels do not occur. The transitions indicated by the 
dashed lines are very weak for large J. 
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for each K value because 

J - K ♦ S, K ♦ S - 1, . . ., |K - S|. 

Since the coupling constant (A) is negative for both electronic states 
involved in the Swan transition, the triplet terms are inverted - that 
is, terms of the triplet components lie in the inverse order of their 
quantum number, {2, which represents the resultant of the electronic 
angular momentum about the intemuclear axis. In a homonuclear molecule 
such as C 2 which has zero nuclear spin, the antisymsetric levels do not 
occur and consequently one component of the A- type doubling does not 
appear in each electronic state. These are the negative levels of the 
state and the positive levels of the 3 ir u state. 

From the selection rules given in table 3.1, P(AJ « -1), Q(AJ ■ 0), 
and R(A J = +1) branches as well as P-form, Q-form, and R-form branches 
are all theoretically possible. Sample transitions corresponding to 
these branches for the K" * 12 level are displayed in figure 3.4. This 
energy level diagram is not drawn to scale and the multiplet splitting 
is greatly exaggerated. By examining the HBnl-London factors which are 
given in table 3.2, one can see that the P-, Q-, and R-form branches 
and the Q branch are extremely weak for large values of J. There- 
fore the transitions corresponding to these branches are given by dashed 
lines in figure 3.4. A close examination of the spectra resulting from 
the transitions on figure 3.4 would reveal that the lines from even- 
numbered K values are displayed to one side of a mean position and 
odd-numbered K values to the other side. This effect in the branches 
is called staggering. 
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Table 3.2 Honl-London Factors for the Swan (d 3 »g - a 3 ir u ) and Fox- 
Herzberg (e 3 irg - a 3 tt u ) Transitions (Schadee, 1964) . These H8nl 
London Factors are Normalized to (2S' ♦ 1)(2J' ♦ 1). 
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3.4.3 The Fox-Herzberg System Transition, e Hg - a 3 ir u 

The Fox-Herzberg system (Fox and Herzberg, 1937) results from 

transitions between the e 3 ir and a 3 tr electronic states. The 

g u 

electronic term energies and vibrational constants required to calculate 
the vibrational energies (equation (3.11)) are given in table 2.1. This 
system occurs in the spectral region between approximately 2300 and 3500 t. 

Figure 3.5 is a Deslandres diagram of Franck-Condon factors of this 
system normalized to the maximum value. The primary Condon locus for 
this system has broadened and its apex has moved diagonally away from 
the (0,0) transition. This feature is typical of band systems in which 
the difference between the equilibrium intemuclear distance (Ar e ) of the 
electronic states is large (see table 2.1). A number of subsidary loci 
can also be seen in figure 3.5. 

The rotational structure of the Fox-Herzberg system is almost 
identical in character to that of the Swan system. Both systems make a 
rapid transition from Hund's case (a) to case (b) as J increases. The 
primary difference between these two systems apart from wavelength, how- 
ever, is that the coupling constant. A, is very small and positive *or 
the upper electronic state of the Fox-Herzberg system. As a result of 
this, the triplet term components of the upper state lie in a normal 
order, that is, terms of the triplet components lie in the order of 
increasing quantum number, n, as shown in figure 3.6. The energy level 
equations, selection rules, and Hbnl -London factors are the same as 
those given previously for the Swan system. Once again the very weak Q 
branch and P-, Q-, and R-form branches are represented by dashed lines. 
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Figure 3.6.- Energy level diagram for the rotational structure of the 
Fox-Herzberg (e 3 ^g - a 3 ir u ) transition, K" * i2. Both electronic 
states belong to Hund's case (b) and the lower state (a 3 ir u ) is 
inverted. The transitions indicated by dashed lines are very weak 
for large J. 
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The staggering effect in the branches can also be observed in the 
spectrum of this system. 

3.4.4 The Ballik-Ramsay System Transition, b 3 Zg - a 3 ir u 

The final triplet system of interest in this work is the Ballik- 
Ramsay system (Ballik and Ramsay, 1959). It occurs when transitions 
take place between the b 3 Z~ and a 3 ir u electronic states. Electronic 
and vibrational term energies are given in table 2.1 and equation (3.11), 
respectively. The approximate spectral range (appendix B) is from 
7100 to 25,000 A. 

Figure 3.7 is a Deslandres diagram of normalized Franck-Condon 
factors for this system. The Condon parabola is narrow as one would 
expect since Ar e for this system is small. 

The lower electronic state is the same as that of the Swan and 
Fox-Herzberg systems and therefore the rotational energy levels are 
given by equations (3.42). For the upper state, b 3 E”, the rotational 
levels are given by 

Fj = B y K(K + 1) + (2K + 3)B y - A - /(2K + 3) 2 B y 2 + A 2 - 2AB y + y(K + l),' 
F 2 = B y K(K + 1), 

F 3 = B y K(K + 1) - (2K - l)B y - A + /(2K + l) 2 B y 2 + A 2 - 2AB y - yK, 

(3.43) 

where Fj, F 2 , and F 3 refer to the levels with J = K + 1, K, and K - 1, 
respectively, and A and y are splitting constants. Again the 3 ir 
lower state is assumed to conform to Hund's case (b) conditions. The 
branches of the 3 Eg - 3 ir u electronic transitions are shown on figure 
3.8. The dashed transitions are extremely weak as one can readily verify 
by consulting the Hbnl -London factors in table 3.3. 
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Figure 3.8.- Energy level diagram for the rotational structure of the 
Ballik-Ramsay (b 3 E» - a 3 ir u ) transition, K M ■ 7. The lower elec- 
tronic state obeys Hund's case (b) conditions and is inverted. 

The P form and R form branches (represented by dashed lines) 
are very weak. 




-54- 


Table 3.3 Httnl -London Factors for the Ballik-Ransay (b 3 E" - a 3 ir ) 
Transition (Schadee, 1964). These Factors are Normalized to 
(2S * + 1)(2J* + 1). 
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3.4.5 The Mul liken System Transition, D l E* - X *E* 

The sisplest of the singlet transitions is the D *E* - X l E* 
transition which is the Mulliken (Mul liken, 1930) system. The term 
energies for the electronic states involved in this transition are 
given in table 2.1, and the vibrational term energies can be calculated 
from equation (3.11) using the vibrational constants in table 2.1. This 

O 

system has an approximate theoretical spectral range of 2100 to 2500 A. 
However, the observed spectral range (see the following paragraph) is 
considerably smaller. 

Figure 3.9 is a Deslandres diagram of the normalized Franck-Condon 

factors for this system. The Condon parabola has degenerated into a 

straight line lying along the principal diagonal of the diagram. This 

phenomenon occurs when the equilibrium intemuclear distances for the 

electronic states are nearly the same. 

Since the selection rules are identical for Hund's case (a) and 

case (b) for 1 E + - *E + transitions, it is immaterial to which case the 
o g 

states belong. The only allowed transitions are those with AJ * AK = ±1. 
Consequently, the band structure which is shown in figure 3.10 consists 
of a single P and a single R branch. The energy levels J* = 0, 2, 

4, 6 and J" ■ 1, 3, 5 missing from figure 3.10 correspond to the anti- 
symmetric levels which as mentioned above do not occur in homonuclear 
molecules with zero nuclear spin. 

The rotational term values are: 

F(J) = B y J(J + 1) - D y J 2 (J ♦ l) 2 ♦ . . ., (3.44) 

and the HBnl- London factors are: 
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Figure 3.9.- Oeslandrcs diagram of normalized Franck -Condon factors 
for the Mull ikon system. Only normalized values exceeding 0.1 are 
included. 
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Figure 3.10.- Energy level diagram for the rotational structure of the 
Mulliken (D - X transition. The enerpr levels J' * 0, 2, 
, . and J" » if 3, 5 correspond to antisymmetric levels and 
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therefore do not occur. 
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It is interesting to note that the Mulliken system is "headless" and 
in this sense is unique among the C 2 systems investigated in this study. 
The appearance of this headless structure is readily understood by not- 
ing the relation between the band origin, v 0 , and the band head, v m , 
spectral locations. 


v - v = - 


<»„. * v .) 2 .. 

> C " • 


(3.45) 


In this equation, B v , and B^,, are the rotational constants for a 
vibrational level v of the upper and lower states, respectively. For 
the Mulliken system, the quantity (B v , - B v „) is very small and conse- 
quently the band head lies a great distance from the band origin and 
corresponds to a large J value. Since the intensity of large J lines 
is approximately zero, the bandhead in these cases is normally too 
weak to be observed. 

3.4.6 The Deslandres-d'Azambuja System Transition, C - A ^ 

The Deslandres-d’Azambuja system (Deslandres and d’Azambuja, 1905) 

results from the C l ar - A V transition. The electronic term 

g u 

energies are given in table 2.1 and the vibrational term energies may be 
obtained by using equation (3.11). This system radiates in the spectral 
region between approximately 2950 and 5400 A (appendix B). 

Figure 3.11 is a Deslandres diagram of the normalized Franck-Condon 
factors for this system. The Condon parabola is narrow which is indicative 
of a small Ar e . 
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Since both electronic states are singlets, it is not necessary 
to distinguish between Hund's cases. From the selection rules 
AJ ■ 0, ±1, it is clear that P, Q, and R branches can all exist. 

Saaple transitions are shown in figure 3.12. The A-type doubling com- 
ponent corresponding to the antisyranetric levels does not appear. In 
A-type doubling as J changes, the symmetric levels are alternately 
above and below the antisymmetric levels. This feature which is greatly 
exaggerated in figure 3.12 produces the "staggering'’ effect in the 
observed spectrum of this system. 

The rotational levels for each state are given by 
F(J) * B y J(J + 1) (A - B y )A 2 - D y J 2 (J + l) 2 ♦ . . . ♦ fj.CJ) (3.46) 
where 4 + (J) is a very small correction to account for the A-type 
doubling. Inspection of the HtJnl-London factors for this transition 
reveals that the Q branch is considerably weaker than the P and R 
branches and its intensity decreases rapidly as J increases. 
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3.4.7 The Phillips System Transition, A - X 1 E + 

An A l * u - X l E + transition produces the Phillips system (Phillips, 
S 

1948). This system has the same lower electronic state as the Mulliken 
system and the same upper electronic state as the Deslandres-d'Azambuja 
system and therefore the electronic, vibrational, and rotational energy 
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Figure 3.12.- Energy level diagram for the rotational structure of the 
Oeslandres-d 'Azambuja (C - A transition. The Q branch 
(dashed line) is considerably weaker than the P and R branches. 
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levels for the lower and upper states may be obtained from sections 
3.4.5 and 3.4.6, respectively. The spectral range of the band origins 
of this system is 5000 to 23,400 A (appendix B). 

A Deslandres diagram of the normalized Franck-Condon factors for 
this system is plotted in figure 3.13. Hie principal Condon parabola is 
rather broad which is typical of a band system with a medium Ar e . Note 
that the apex of the principal parabola is still located at (0,0). A 
subsidiary parabola also appears in the diagram. Its presence is also 
typical of medium (and large) Ar e 's. 

For this electronic transition, the distinction between Hund's case 
(a) and case (b) is again needless. The selection rules AJ ■ 0, ±1 
limit the transitions so that only simple P, Q, and R branches occur. 
These transitions are shown on the energy level diagram in figure 3.14. 

The odd J" levels are antisymmetric and do not appear. 

The H5nl- London factors are: 

J f A 1 

Branch ^T'A" 

P (1/2) (J - 1) 

Q (1/2) (2J ♦ 1) 

R (1/2) (J ♦ 2) 

3.4.8 The Freymark System Transition, E l E* - A 1 » u 

The Freymark system (Freymark, 1951) results from transitions between 

the E and A 1 ir„ electronic states. The electronic term energies 
g u 

are given in table 2.1, and the vibrational term energies can be obtained 
from equation (3.11). The band origins (appendix B) for this system 
range from 1950 to 25S0 A. 
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Figure 3.13.- Deslandres diagram of normalized Franck-Condon factors 
for the Phillips system. Only normalized values exceeding 0.1 are 
included. 
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Figure 3.15 is a Deslandres diagram of the normalized Franck- 
Condon factors for this system. The Condon parabola is narrow as expected 
since the Ar e for this band system is small. 

The rotational energy levels for the upper and lower states can be 
obtained from equations (3.44) and (3.46), respectively. Sample transi- 
tions are shown on the energy level diagram in figure 3.16. 

The HBnl -London factors for this system are: 

ji a* 

Branch 

P (1/2) (J + 1) 

Q (1/2) (2J + 1) 

R (1/2)J 

3.5 Shock-l\ibe and Thermochemistry Theory 

This section reviews the parts of shock-tube and thermochemistry 
theory necessary to calculate the equilibrium concentration, ^ uv » ^ » j < • 

This quantity is an input parameter for equation (3.1) and must be 
accurately known in order to calculate the rotational line intensities 
which are used in the data analysis technique. 

3.5.1 Thermally Perfect Gases in Equilibrium 

The thermodynamic properties of most gases at moderate temperatures 
and pressures are well approximated by an equation of state of the form: 

PiV = NjkT. (3.47) 

By definition, a thermally perfect gas is one which obeys this equation 
of state. The gas mixture in this study is assumed to consist of 
chemically reacting, thermally perfect species. 
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9 0 . 50 - 0.74 
9 0 . 25 - 0.49 
9 0 . 10 - 0.24 


Figure 3.13.- Deslandres diagram of normalized Franck-Condon factors 
for the Freymark system. Only normalized values exceeding 0.1 are 
included. 
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e , r + 

K'-J* E 



Figure 3.16.- Energy level diagram for the rotational structure of the 
Freymark (E - A *ir u ) transition. Odd J' levels are antisym- 
metric and do not occur. 
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The basic condition for theroochemical equilibrium is that at a 

constant temperature and pressure the concentrations of the various 

reacting species must be such as to minimize the Gibbs free energy of 

the gas mixture. It should be noted that conservation of mass and 

charge place certain constraints on this minimization. Gibbs free 

energy, F, for a single species thermally perfect gas is given by 

Vincenti and Kruger (1965) as 

F - NkT(tn N - An Q) (3.48) 

where N and Q are the particle concentration and partition function, 

respectively. Therefore the total Gibbs free energy of a gas mixture 

consisting of j species is 

F = kT 5^ N.[*n N. - in(Q).]. (3.49) 

y J J J 

Applying the free-energy minimization principle to this relation while 
the temperature and pressure are held constant gives 

dF = kT 23 dN,[in N ■ - *n(Q).] ♦ kT 23 N * d(tn N.) * 0. 

4 J J J i J 3 

3 3 (3.50) 

Since d(An Nj ) = dNj/Nj , the last summation in this equation vanishes 

and the final condition for equilibrium is: 

23 dN,[m N. - tn(Q) . ) » 0 (3.51) 

j 3 3 J 

for all values of dNj which satisfy the mass and charge conservation 
constraints. Applying equation (3.51) to the general chemical reaction, 

“1*1 ♦ “2^2 ♦ • • • + eyc^ + °£ + i x t+1 ♦ • • • ♦ °j x j* (3.52) 
one obtains 

Si 5 ® - n (Q,)°^, (3.53) 

v i 3 j 3 
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from which the equilibrium species concentrations can be calculated 
provided the partition functions are known. In the above equations, 
the stoichiometric coefficients, o^, are negative for reactants and 
positive for products. Following the usual convention, the equilibrium 
constant is defined by equation (3.53). 

The partition function Q is defined by the expression 

(3.54) 


— ‘ E «/ kT 
q ■ Ev 

n « 


where the n summation extends over all known energy levels. Since 

molecules (and atoms) possess various modes of energy, the total energy, 

E, of a molecule can be closely approximated by 

E = + E +E..+E, . (3.55) 

trans rot vib elec 

By substituting this expression into equation (3.54), one finds that the 
total partition function Q is expressible as the product of individual 
partition functions for each mode of energy, that is, 

^ * ^trans * *Vot * %ib * ^elec * ^trans * ^int* (3.56a) 

In order to calculate the thermochemical properties of a gas mixture, 

it is necessary that each species be referenced to a common zero of 

energy. For this requirement, the partition function becomes 

-h^kT 


^ * ^trans * ^int * e 


(3.56b) 

-trans xnt 

where hj° is the heat of formation of the species i and represents 
the energy difference between the lowest quantum state of the species i 
and the zero of energy. 

The translation partition function is given by (Herzberg, 1950) 

V3/2 


trans 


, / 2xmkT V 

V » 2 / 


(3.57) 
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where m is the mass of the particle and V is the volume of the 
system under consideration. 

Hie energy Eji of an individual molecular energy level in an 
electronic state n is given by (section 3.1) 

En - hc[T e + G(v) ♦ F(J)] n (3.58) 

where T e is the term energy of the electronic state referenced to the 
lowest level (v » 0, J * 0) of the ground state; G(v) is the vibrational 
energy of the vth level above the value T e ; and F(J) is the rota- 
tional energy of the Jth level above the value T e ♦ G(v). With this 
expression, the internal partition function for diatomic molecules is 


r.yT (2 ' S o,A )( 2S * 1)(2J * l) 

^ ' vrr ’ ° 


exp 


{(' w) lT e + G(v) . + 

(3.59) 


where (2-6 .)(2S + 1) (2J + 1) is the total statistical weight of each 

0, A 

level, a is a symmetry factor which approximates the nuclear spin 
effects of homonuclear molecules - it is set equal to 2.0 for homonuclear 
molecules and 1.0 for all others. 

For atoms, since vibrational and rotational energy modes are absent, 
the total partition function becomes 

o “ WW < 3 - 6 °> 

where Q trans is again given by equation (3.57) and Q e j ec * Qj nt is 
given by 

Qint * £ *n exp(-En/ k T). (3.61) 

n 

In this equation, g n is the statistical weight of the level n. 



-71- 


The equilibrium concentrations, N uv » j » a * * can now be calculated 
using equations (3. S3) and (3.57) with equations (3.S9) and (3.61) for 
molecular and atoms, respectively. Before this can be accomplished, 
however, the thermodynamic state of the gas must be obtained. This 
procedure is outlined in the following section. 

3.5.2 Shock -Tube Theory 

The gasdynamics of shock waves is well documented in many sources. 
Among these are the textbooks by Liepmann and Rosko (1957) and Vincenti 
and Kruger (1965). The following brief discussion was taken from these 
references . 

For a flow which contains a normal shock wave (see fig. 3.17) the 
conservation equations of fluid flow written in terms of shock-fixed 


coordinates are: 



Mass: p^u^ * 

p 2 u 2 

(3.62) 

Momentum: P t 

+ PlUj 2 « P 2 ♦ p 2 u 2 2 

(3.63) 

Energy: hi + 

(l/2) Ul 2 = h 2 + (l/2)u 2 2 

(3.64) 


Shock-heated 

Initial 

test gas 

conditions 

© 

__ © 


1 


SHOCK 

FRONT 

Figure 3.17.- Sketch of normal shock flow in a constant area shock tube. 
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where the subscripts 1 and 2 denote conditions upstream and downstream 
of the shock, respectively. Equations (3.62) -(3.64) assume a constant 
area shock tube. The shock velocity is equal to ui. These equations 
are equivalent to the Rankine-Hugoniot equations and give the jump 
conditions that exist when a shock wave is represented by a discontinuity. 
If the initial conditions 1 are known, the set of equations (3.62)- 
(3.64) can be solved by assuming the gas to be in thermodynamic equilib- 
rium on each side of the shock wave. This assumption adds the state 
equation 

p = p(P, h) (3.65) 

and one then has four equations and four unknowns. Since the relation- 
ship in equation (3.65) is normally not given in analytical form, the 
usual procedure is to solve the set of normal shock relations using an 
iterative technique described by Vincenti and Kruger. In this technique, 
equations (3.63) and (3.64) are rewritten in the form: 

P2 * Pi + Piui 2 [l - (p 1 / P 2 ) ] » (3.66) 

h 2 - h! ♦ (ui 2 /2)[1 - (Pi/p 2 ) 2 ). (3.67) 

Using these equations and an estimated value of (P 1 /P 2 ). trial values 
of P 2 and h 2 are obtained. These results are used in the equation of 
state to obtain a value of P 2 and thus a revised value of (Pi/P 2 )> 

This procedure is repeated until the initial and final values of (P 1 /P 2 ) 
are identical. With P 2 , h 2 , and P 2 known, all of the other thermody- 
namic parameters as well as the equilibrium concentrations N uv , a ,j, 
can be calculated. 

The application of the equations, reviewed in this section and 
the preceding one, to the present study is presented in the next section. 
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3.5.3 Thermochemistry Calculations 

The computer program used to calculate the equilibrium thermo- 
chemical properties of the shock-heated gas mixtures was developed by 
Horton and Menard (1969) . In this program, the thermodynamic properties 
of the shocked gas are determined by an iterative procedure similar to 
that described in section 3.5.2. A final temperature T 2 is specified 
and (pi/p 2 ) is varied until the Rankine-Hugoniot equations are 
satisfied. The solution is assumed to be converged after n iterations 
if 

(Pl/P2) n+1 - (Pl/P2) n | 


(Pl/P2) 


n 


< 10 


-3 


(3.68) 


With the thermodynamic properties completely determined, the equilibrium 


species concentrations can then be calculated from equation (3.53), 

f/2nm.kT\ 3 / 2 


0 r/2nmkT\ 3 / 2 

k n ■ " <v * (Ql) jJ ’ 


(3.69) 


and the conservation of mass equation, 

"j ’ ?Vi 


(3.70) 


where is the number of atoms of the jth numbered element in the 

molecule of the ith numbered species. 

In this work, the internal partition function of molecules 


Qint ?7T 0 


«p{(- lifK * otv) * F(j)i l„ 

(3.71) 


was approximated as described below. 

The rotational partition function was approximated (see Herzberg, 
1950) by 
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£ (2J ♦ l)exp[- 



(2J ♦ 



B y J(J ♦ 



(3.72) 


kT I 

= hcB y ) 

where B = B - a (v + (1/2)]. 
v e e 

The vibrational energy within each electronic state was assumed to 
be given by 

G(v) = w [v + (1/2)] - u x [v ♦ (1/2) ] 2 ♦ u V [v ♦ Cl/2)] 3 ♦ « zjv ♦ (1/2)]**. 

c c c c “ C " 

(3.73) 


Therefore, the internal partition function for diatomic molecules was 
calculated from 


^int 


?( 


( 2 - « 0 , A H2S + 1) 


' ' v *n 


(3.74) 

The summation is over all electronic levels of importance and over all 
vibrational levels in each electronic state until 


e xp.L:. hc G(v.VhII..- < 10 -3. (3.7 

2 exp[-hcG(v)/kT] 
v 

Appendix D contains a list of the species included in the thermo- 
chemical calculations of mixtures of argon and acetylene. The energy 
levels and other significant constants pertinent to the calculations 


are also included. 
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3.6 Synthetic Spectrum Data Analysis Technique 

In this section, the data analysis technique which incorporates 
synthetic spectrum calculations based on equation (3.1) is described. 

The assumptions and approximations to the theory presented earlier on 
which these calculations are based are discussed. Finally, the details 
of the reduction of absolute measurements of spectral emission to 
electronic transition moments, E|R e /ea 0 | 2 , are disclosed. 

3.6.1 Synthetic Spectrum Calculations 

Synthetic spectra were calculated using the computer program 
described by Arnold, Whiting, and Lyle (1969) and Whiting, Arnold, and 
Lyle (1969). This program produces a spectrum by accounting for the 
contribution of each rotational line (equation (3.1)) of interest in 
an electronic transition. 

The integrated spectral intensity due to spontaneous emission for 

each rotational line was previously shown to be given by 

16 x 10' 7 _ JtA , 

E t 3 3(2 - 6 a,)(!2S' + 1) (2J' ♦ 1) N uv'A'J» Z l R e (r v V' 5 ^ 2q v , v" S J"A M * 

O j A 

(3.76) 

The input parameters used by the computer program to compute E^ 11 are 
discussed in the following paragraphs. 

Hie central wavelength of a rotational line is calculated from 

v - v/c - v o (v', v") + F* (J') - F H (J"), cm’ 1 (3.77) 

where the band origin v 0 (v', v") is given by 

v (V, v") » (T ' - T ") ♦ (G'(v') - G"(v") J • (3.78) 

O 6 c 

The expressions for F(J) listed in section 3.4 are approximated in the 
program by 
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F(J) - B y J(J ♦ 1) - D y J 2 (J ♦ l) 2 . (3.79) 

This expression is used for all singlet transitions and all triplet 
transitions in which spin splitting is negligible with the exception 
of the AA ■ 0 triplet cases. For these transitions, the central com- 
ponent of the triplet is used as an effective line wavelength and F(J) 
is approximated by 

F(J) - B y [J(J ♦ 1) ♦ 4Z 2 ] - D y [J ♦ Cl/2)] 4 (3.80) 

where 

2 . AY(Y . 1) - (A/9) - 2J(J . 1) (Jill) 

3[A 2 Y(Y - 4) . (4/3) ♦ 4J(J ♦ 1)) 
and 

Y » A/B y . (3.82) 

The quantity G(v) is approximated by 

G(v) - Ue [v * (1/2)] - u) e x e [v ♦ (1/2)] 2 + <u e y e [v ♦ (1/2) ] 3 ♦ ^(v + (1/2)]“. 

(3.83) 

Molecular constants such as T , o> . u x , etc. required as program 
inputs were taken from Rosen (1970) and Herzberg (1950) and are listed 
in table 2.1. 

The Franck-Condon factors, q v , v „, used in the calculations were 

discussed in section 3.2 and are listed in appendix C. 

J' A' 

The expressions for the Honl-London factors, listed in 

section 3.4, are approximated in the program as follows: 

1. Parallel transitions (AA « 0): These transitions (Swan, Fox- 

Herzberg, Mulliken, and Deslandres-d'Azambuja systems) are characterized 
by strong P and R branches and very weak Q, P-form, Q-form, and R- 
form branches. The H8nl -London factors for these transitions are 
approximated in the program by: 
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Branch <J J H A U 

P J 

R J ♦ 1 

These strength factors are exact for the * £ - transitions (section 
3.4.5). In this assignment, the weak branches are assumed to have 2 ero 
intensity. However, the intensity actually possessed by these transi- 
tions is incorporated into the P and R branches because the sum of the 
HBnl- London factors for a given J value is held equal to 2J + 1. 

2. Perpendicular transitions (AA = ±1): These transitions 

(Ballik-Ramsay, Phillips, and Freymark systems) are characterized by 
strong P, Q, and R branches and very weak P-form, Q-form, and R-form 
branches. The HOnl- London factors for these transitions used in the 
program are: 

J' A' 

Branch ^"A" 

„ ( J ± A) (J - 1 t A) 

v 23 

(J A A)(J ♦ 1 * A) (2J ♦ 1) 

4 2J(J ♦ 1) 

d (J + 1 l A) (J ♦ 2 l A) 

K 2(J ♦ 1) 


where the upper sign is for AA » +1 transitions and the lower sign is 
for AA = -1 transitions. As discussed above, the intensity of the 
weak branches is accounted for in the P, Q, and R branch intensities.. 

In the program, the number of particles in the upper state is 
calculated from 


N 


uv'A'J' 


N du(2J* + 1) 


exp 1 


•he 

nr 


T • ♦ C ' (v ' ) + F ' (J 


2 1 ) 


(3.84) 
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where du is the electronic multiplicity of the upper state. The 
partition function, Q (see section 3.5.1), is approximated by 


Si 


(2 - « 0>a )(2S + 1) 




(3.85) 


where the summation is carried out over all electronic levels of 
importance for the temperatures considered and over all vibrational 
levels within each electronic state until either 


exp[-hcG(v)/kT] < 1Q - 3 
£ expf-hcG(v)/kT] ' 


v 

or the equation repressing the vibrational energy (equation (3.83)) has 

reached a fictitious peak, that is, G(v + 1) - G(v ). 

r v max ' v max' 

Once the integrated intensity of each line has been calculated from 
equation (3.76), this intensity is then distributed in the spectrum by 
an approximate Voigt profile which is given by 




[(* - Xcl)/W v ]2.2 5 

(3.87) 


The spectral intensity at the line center is 

E * U 

E X * i , (3.88a) 

T. W y [1.065 + 0.447(W L /W V ) + 0.058(W L /W y ) 2 ] 

where 


W y = W l /2 + /Sf L /2) 2 + W g 2 . 


(3.88b) 



-79- 


In the above equations, Wy, W, , and W are the Voigt, Lorentzian, and 

© 

Gaussian profiles at half-height. The Gaussian line profiles at half- 
height were calculated from: 

W = 7.16 x 10 -7 X/f/M , l (3.89a) 

8 

and are related to the Lorentzian profiles at half-height by 

W g = v'tnT/a (W L ) . (3.89b) 

In the above equations, T is the temperature of the radiating gas, M 
is its molecular weight, A is the wavelength of the radiation, and a 
is a parameter which is usually determined from experimental measi.e- 
ments. However, since the measurements were made in the optically thin 
regime (see Chapter 4), the exact line profile used in the calculations 
is not critical and therefore the value of a used is not either. As 
an example, the integrated intensity of the Av = 0 sequence of the Swan 
system was calculated for values of a ranging from 1.0 to 0.001. The 
maximum difference in integrated intensity was only 2.7 percent. 

The contribution from each line is then summed yielding a synthetic 
spectrum. 

3.6.2 Data Reduction Technique 

The data reduction technique used the synthetic spectrum calculations 
discussed in the previous section. A synthetic spectrum was calculated 
by assigning a value of unity to the Z|R e /ea 0 | 2 value. This spectrum 
is then folded with the measured experimental instrument function (fig- 
ure 4.5) to give a predicted signal output. The ratio of measured to 
theoretical signal output yields tKe experimental value of E|R e /ea 0 | 2 . 

Before comparison with the theoretical signal, ali of the 
experimental signals were first corrected to account for the presence of 
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weak oblique shocks in the radiometer viewp&th (section 4.2.4). This 
correction was approximately 2 percent. For the measurements taken in 
the infrared region., the experimental signal was also corrected for 
atomic line radiation occurring in the bandpass of interest. The inte- 
grated intensity due to spontaneous emission of a single atomic line is 
calculated from 

NhvA, 

E » - u - ^ - - u * , W/cm 3 »sr. (3.90) 

All of the allowed transitions for Arl and Cl listed by Wiese et al. 

(1966) were included in the calculations. The Einstein A coefficients, 

A # , for these lines were also taken from Wiese et al. 
ux 

The partition function used in the calculation of the number of 
atoms in the upper state, N u , was identical to the one used in the com- 
putation of the species concentrations. As was discussed for molecules, 
the integrated intensity of individual atomic lines is summed and then 
folded with the instrument functions to give a predicted signal output. 
This output is subtracted from the measured output before comparison 
with the theoretical signal based solely on molecular radiation. The 
largest correction for atomic line radiation was 1.2 percent. 

All theoretical calculations are based on optically thin approxi- 
mations. This assumption was justified by comparing synthetic spectra 
based on the present z|R e /ea 0 | 2 measurements with black-body curves. 

The black-body intensities were always at least an order of magnitude 
greater than the peak spectral intensities. Line widths used in these 
calculations were derived from the formulas given by Arnold et al. 

(1969) and neglected instrument broadening. 
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Chapter 4 

Description of Experimental Technique 

This chapter contains a detailed description of the experimental 
facility and instrumentation used in the present work. The first 
section describes the shock-tube facility and its operational procedures. 
Section 4.2 covers the instrumentation and calibration procedures used 
in the experimental program to obtain reliable measurements. The third 
section (4.3) outlines the initial test conditions, and section 4.4 
presents samples of the experimental data. Survey spectra and oscillo- 
grams of the spectral emission for each of the band systems measured are 
shown, 

4.1 Shock-Tube Facility and Operation 
4.1.1 Idealized Shock-Tube Flow 

Several advantages are gained by using a shock tube to provide a 
heated sample of gas. The gas sample is raised nearly instantaneously 
to a high temperature, pressure, and density; the gas equilibrates 
rapidly and maintains this state for a short period of time; and the 
state of the shocked gas can be calculated very accurately. A brief 
description of the wave pattern and state of the gas inside a conven- 
tional constant-area shock tube is given below. For additional details 
of this subject, the reader is referred to the work of Gaydon and 
Hurle (1963) and Zel'dovich and Raizer (1966). 

A sketch of a simple shock tube and an accompanying distance- 
time (x-t) diagram showing the paths of the waves is shown in figure 
4.1. Initially there are two regions of gas at vastly different 
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Figure 4.1.- A distance-time diagram showing the progress of the 
flow in a conventional shock tube. 

pressures separated by a diaphragr. The gas in the low pressure side 
(conditions "1") is the experimental sample, while the gas in the high 
pressure side (conditions "4") is known as the driver gas. 

Upon bursting of the diaphragm which is assumed to open 
instantaneously, the driver gas expands into, the test gas, compressing 
it and forming a shock wave which propagates through it. Because of 
this compression, the test gas behind the incident shock wave (conditions 
"2") has undergone a substantial increase in tempera cure, pressure, and 
density. Since the initial conditions, "1", and the composition of the 
test gas are known, the velocity of the shock wave is all that is 
required in order to calculate the conditions "2" immediately behind the 
shock. This numerical procedure uses the Rankine-Hugoniot equations 
and was outlined in section 3.5. 
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As the shock wave propagates through the test gas, an 
expansion or rarefaction wave simultaneously travels into the high 
pressure section at a velocity equal to the local speed of sound. The 
surface separating the two gases is known as the contact surface. 

The above description is an idealized picture of shock-tube flow. 

In reality, the shock tube performs somewhat differently, for example, 
the contact surface is usually not an abrupt boundary since diffusion 
and other processes tend to mix the particles from both regions. This 
mixing in addition to the formation of a cool boundary layer on the 
shock-tube wall behind the incident shock tend to reduce the test time, 
that is, the length of the test gas free of contamination divided by 
the shock velocity. Ideally the shock wave propagates at a constant 
velocity while in normal operation most shock waves tend to slow down 
due to viscous effects as the shock moves along the tube. 

4.1.2 Shock-Tube Facility 

The facility used in the experimental program is a combustion- 
driven shock tube and is located at Ames Research Center, Moffett Field, 
California. A photograph showing an overall view of the facility is 
presented as figure 4.2. The driver section (combustion chamber), 
which is the largest diameter portion of the shock tube, is shown in 
the foreground a p 1 is 1.27 m long. It has an inside diameter of 16.3 cm. 
A preformed stainless-steel diaphragm (known as the primary diaphragm 
and not shown in the photograph) separates the driver and driven tubes. 
The driven tube is 10.3 m long and has an inside diameter of 16.8 cm. 

This tube, shown in figure 4.2 extending from the center of the photo- 
graph to the upper left-hand corner, is coupled to a stainless-steel 




Figure 4 . 2 .- Photograph of the combustion-driven shock-tube facility located at Ames Research Center, 

Moffett Field, California. 
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dump tank (extreme left of photograph) which is approximately 1 m in 
diameter and 2 m long. A flat thin prescored stainless-steel 
diaphragm separates the dump tank from the driven tube. The driven 
tube contains a test section whose centerline is 9.2 m from the primary 
diaphragm. The test section consists of four ports 5 cm in diameter 
and 90° apart. These ports can be used for obtaining either photographic, 
spectrographic, or radiometric data. Six additional ports at various 
locations along the shock tube and nominally 1 m apart are used to 
accommodate positive ion sensors which are used in determining the 
velocity of the shock wave. 

The shock tube represents a relatively tight vacuum system A 
Welch two-stage rotary vacuum pump is used to routinely evacuate the 
driven tube to a pressure of less than 10" 3 torr. The leak rate of 
this section is nominally 5 * 10" 5 torr/min. 

A continuous flow-through system for loading the test gas has 
been incorporated into this facility. The initial pressure in the 
driven tube is set by adjusting the bleed-in and pump-out rates until 
the predetermined pressure is obtained. Pure test gas flows through 
the tube until approximately 6 sec before the shock tube is fired when 
two solenoid operated valves at the bleeJ-in and pump-out locations 
terminate the process. Mock tests in which the shock tube was not 
fired verified that the gas pressure in the driven tube did not change 
as a result of terminating the gas flow. In this work, both the 
flow-through loading technique and simply loading the tube to the 
desired pressure were used. The radiometric and spectrographic results 
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fron these tests were essentially identical. Therefore to minimize 
costs, most of the tests were performed by loading the test gas to a 
desired pressure and firing. 

The energy for initiating combustion is provided by a 20 MFD 
capacitor charged to S KV. In order to promote uniform combustion in 
the combustible gases, this energy is discharged into a 10 mil tungsten 
glow wire positioned at the centerline and extending the entire length 
of the combustion chamber. 

4.1.3 Shock-Tube Operation 

In the following paragraphs, the procedures for cycling the shock 
tube through a typical test are presented. A thorough cleaning of the 
entire shock tube with acetone wetted cloths was the first step. All 
port windows including those in the splitter plates (see next section) 
were either cleaned or replaced. At this point the radiometric and/or 
spectrographic instrumentation was carefully aligned about the test 
section. Wavelength calibrations and band-pass measurements of the 
radiometers were also performed. 

Both the glow wire and a powder-driven squib-ignited punch which 
is used to open the primary diaphragm were installed. The shock tube 
was then assembled using diaphragms which divide the tube into three 
sections: combustion chamber, driven tube, and dump tank. All three 

chambers were then evacuated. The driven tube was then purged by fill- 
ing it with test gas to a pressure which was nominally twice the initial 
pressure, Pj, and reevacuating. It should be noted that a dry-ice 
cold trap was used on the test gas bleed-in line to remove moisture. 

The test gas was then loaded to the desired initial pressure. 
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Next, the combustion chamber was loaded with a combustible 1:3:8 
mix of oxygen, hydrogen, and helium. With the completion of this task, 
the shock tube was ready to be fired. This was accomplished by dis- 
charging the energy from the 20 MFD capacitor into the glow wire. This 
process ignited the combustible gases which heated the mixture and 
caused the pressure to increase approximately six fold. Shortly after 
the peak pressure was achieved in the combustion chamber, the punch 
pierced the primary diaphragm at its center and the weakened diaphragm 
was forced open by the high pressure gas. This gas rushed into the 
driven tube and created a shock wave which propagated through the gas 
heating it. Ion sensors spaced along the tube detected the arrival of 
the shock and recorded information from which the shock velocity could 
be deduced. The shock wave passed through the test section where the 
radiation from the shock-heated gas was measured and then reflected off 
of the secondary diaphragm. The initial contact weakened the diaphragm 
and the subsequent higher pressure opened the diaphragm and the gases 
flowed into the dump tank ending the test. A description of the instru- 
mentation used to determine both the thermodynamic and spectroscopic 
properties of the gas mixture is given in section 4.2. 

4 . 2 I nstrumentation 
4.2.1 Pressure Measuring Devices 

A 0-5000 yHg Stokes McLeod vacuum gage was used to monitor the 
minimum driver tube pressure and leak rate. A 0-20 torr Wallace and 
Tieman vacuum gage was used to set the initial pressure in the driven 
tube. Both gages were calibrated by direct comparison with a reference 


standard traceable to the National Bureau of Standards. A 
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piezo-electric pressure transducer mounted in the breech plug of the 
combustion chamber was used to measure the driver pressures. 

Figure 4.3 is an oscillogram showing the pressure history in the 
combustion chamber. The oscillogram sweep was started at time-zero. 
Approximately 35 msec later (indicated by the first small pulse on the 
pressure record) , the energy was discharged into the glow wire. The 
combustion process heated the gas and caused a rapid increase in pres- 
sure in the chamber. Some 65 msec after time- zero and after the pres- 
sure had reached an essentially constant value, the punch was fired 
(indicated by the second small pulse on the oscillogram). The sudden 
decrease in pressure is caused by the arrival of the expansion wave at 
the breech end of the combustion chamber. 

4.2.2 Shock Detectors 

Shock speeds were determined from time-of-arrival measurements made 
at six locations along the driven tube by positive ion sensors. The 
output of each detector was connected to a Hewlett-Packard 100 mHz elec- 
tronic chronograph to record the shock time of arrival . These detectors 
were quite reliable and allowed the determination of shock speed to an 
accuracy of about 0.10 percent. A plot of the variation of shock speed 
with distance for a typical test is shown in figure 4.4. 

4.2.3 Spectrographic Instrumentation 

Spectrographs were used in some of the tests to obtain survey 
spectra of the shock-heated gas for identification purposes. The 
spectrograph viewpath was normal to the shock-tube axis. Three 
spectrographs were used altogether. Only one was used in any particular 
test and all of them employed an opaque shutter, similar to the one 
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Figure 4.3.- Oscillogram of the pressure history in the combustion chamber. 
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Figure 4.4.- Shock speed variation with distance for a typical test. 
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described by Borucki et al. (1969), to exclude the radiation from 
the driver gas. This shutter was explosively driven and had a closing 
time of 35 ysec. 

To survey the ultraviolet portion of the spectrum, a Huet spectro- 
graph model UV 24 was used. This instrument is a dual quartz prism 
spectrograph which has a f/3.5 aperture ratio and an 0.24-m focal 
length. The linear dispersion varies from 7.S A/mm at 2150 A to 
100.0 A/ mm at 4500 A. 

The second spectrograph, used to cover the ultraviolet and visible 
portions of the spectrum, was a model 75-000 Jarrell-Ash grating spec- 
trograph. This instrument has a f/5.6 aperture ratio, a 0.75-m focal 
length, and a dispersion of approximately 20 A/mm. The 600 lines/mm 
grating is blazed for 5000 A. 

The final spectrograph, used in the infrared, was a stigmatic 

O 

f/2.5 grating instrument with a dispersion of about 50 A/mm. The grat- 

O 

ing has 600 lines/mm and is blazed for 7S00 A. This instrument was 
designed by Borucki (1970). 

Sample spectra of the shock -heated gas obtained with these 
instruments are presented in section 4.4. 

4.2.4 Radiometric Instrumentation 
A. General description of radiometers 

Absolute intensities of selected spectral features of the radia- 
tion were measured behind the incident shock wave as a function of time 
by absolutely and spectrally calibrated narrow band-pass radiometers. 
The radiometers consisted of photomultiplier tube detectors coupled to 
Bausch and Lomb 1/4-m monochromators or to interference filters which 
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acted as dispersive elements. The interference filters (used in the 
infrared) incorporated a Corning 2-64 sharp cut-off filter to exclude 
light below 6400 X. Below 3000 A, an additional narrow band-pass 
interference filter was used in connection with the monochromators to 
exclude stray light. 

Figure 4.5 is a sketch of the radiometric instrumentation arrange- 
ment about the shock-tube test section. As illustrated in the figure, 
the radiometers viewed the flow through narrow (0.5 mm) slits in 
splitter plates. Quartz windows were flush mounted on the surface of 
the splitter plates directly in front of the slits so that the entire 
surface of the splitter plate was completely smooth. The slits were 
normally positioned so that the portion of the splitter plate mounted 
on the opposite wall visible to the radiometer was blackened. Addi- 
tional radiometer units consisting of one or five radiometers and 
mounted on the top and opposite side of the shock tube are not shown 
in figure 4.5 for clarity. Thus a total of four splitter plates were 
mounted in the shock tube with the distance between the plates maintained 
at 11 cm. The actual number of radiometer units used each test varied 
depending on the strength of the band system being measured. Measure- 
ments of the Freymark and Mul liken systems, for example, were made with 
single close-coupled radiometers mounted on the side and top of the tube. 

The splitter plates were used to exclude from the radiometer 
viewpath the relatively thick variable temperature boundary layer on the 
shock-tube wall, replacing it by the thinner splitter-plate boundary 
layer. These plates and their accompanying struts were designed accord- 
ing to the following criteria: (1) the strut had to be long enough 
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Figure 4.5.- Sketch of a two-channel radiometer unit coupled to the shock 
tube. The radiometer bandpass featured in the inset is typical of 
those resulting from the 1/4-m mono chroma tor dispersive elements. 
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to completely reject the wall boundary layer, and (2) the splitter 
plate had to be large enough to keep the flow disturbance created by 
the strut between the shock-tube wall and the splitter plate. Photo- 
graphs of the flow over the splitter plate assembly (shown below) 
verified that the splitter plates succeeded in meeting these objectives. 

While the function of the splitter plates was to exclude the wall 
boundary layer from the radiometer viewpath, it is well known 
(Gruszczynski, 1968 and Arnold and Nicholls, 1970) that splitter plates 
cause weak oblique shocks to form on top of the plates which appear in 
the viewpath. Figure 4.6 is a self-luminous photograph of the flow over 
one of the splitter plates and a diagram illustrating the flow processes. 
This photograph was taken with an Abtronics image converter camera 
(model NR-2HS) . The exposure time is 0.5 psec. With careful examina- 
tion of the photograph, an oblique shock originating at the leading 
edge of the splitter plate can be seen. It is also apparent that there 
is no flow over the sides of the splitter plate. 

Equilibrium calculations, using the oblique shock relations and 
the angle of the shock measured from photographs similar to figure 4.6, 
were performed to determine the effect of the presence of these shocks. 
The oblique shock relations can be obtained from the normal shock rela- 
tions in section 3.5.2 by simply replacing the velocity uj by the 
quantity Uj sin 9 which is the component of velocity normal to the 
oblique shock whose angle is given by 0. These calculations revealed 
that the oblique shocks on the splitter plates increased the intensity 
observed by the radiometers by only 2 percent. 



NORMAL SHOCK WAKE 



Figure 4.6.- Sketch and self-luminous photograph of the flow across a 
splitter plate. The strong shock below the splitter plate is 
clearly visible while the weak shock above the plate (depicted in 
the sketch) is barely visible. 
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B. Calibration procedures 

Craig and Oavy (1969) have shown that the radiometer output V(A) 
may be expressed by the following equation: 


V(A) = G(A) 


r 


T(A)H(A)dA, 


(4.1) 


where G(A) is the sensitivity of the system or the calibration constant, 

T(A) is the normalized spectral response of the radiometer, and H(A) 

is the spectral irradiance incident on the radiometer. For the response 

* 

to a calibration source, the quantity H(A) is replaced by the spectral 
irradiance from the calibration source, L(A), and equation (4.1) becomes 


f 


V C (A) = G(A) I T(A)L(A)dA. 


(4.2) 


The ratio of radiometer response from the test to that from the 
calibration source is 


V t (A) 

V C (A) 


r 


T(A)H(A)dA 


r 


(4.3) 


T(A)L(\)dA 


H(A) is the quantity which is to be determined in the present tests. 
For most radiators, it can be written as 


H(A) = H 0 h(A), (4.4) 

where Hq is a constant which gives the absolute strength of the 
radiator and h(A) is the relative spectral irradiance. Since molecular 
bands are the dominant radiators in this work, the quantity Hq can be 
replaced by the electronic transition moment, £|R e /ea 0 | 2 . The relative 
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spectral irradiance, h(x)» was previously described in section 3.6. 
Substituting these quantities into equation (4.3) and rearranging, one 
obtains 


, „ I T(X)L(»dA p T()) - . 

^l ea o| [v c (X)J 


f T(X)h(X)dX 

*o 


v T (x) 


w 

G C (X)J T(X)h(X)dX 


(4.5) 


Therefore, the zlRg/eagl 2 value can be determined from the radiometer 
output during a test if the calibration constant, G C (X), and the spectral 
response, T(X), of the radiometer is known. 

Each radiometer was calibrated absolutely against an NBS tungsten 
ribbon filament lamp standard or a quartz-iodine tungsten-coil filament 
lamp to determine G C (X). The calibration source was used as an infinite 
wall source. The spectral irradiance of the standard lamp was extended 
below 2500 X by measuring the relative output of the lamp between 2200 X 
and 3000 X (at 100 X intervals) and matching the curves at 3000 X. This 
procedure resulted in a difference of only about 5 percent at 2500 X 
between the normalized curve and the calibration value provided with the 
lamp. It should be noted that the optical path, monochromator settings 
(where applicable), and photomultiplier operating voltages were identical 
for both calibration and experiment. 

A typical normalized spectral response, T(X), of a radiometer using 
a 1/4-m Bausch and Lomb monochromator as a dispersing element is shown 
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in the inset on figure 4.5. This response, often referred to as a 

O 

bandpass, was obtained by scanning the 5461 A line of Hg and normaliz- 
ing the output. Both the spectral response, T(A), and the center wave- 
length of the response curve are critical to the understanding of the 
experimental measurements. Consequently, the spectral response, which 

' O 

gives an accurate wavelength calibration (± 1.5 A) was measured before 
every third test. The only response changes detected during this 
procedure were small ones. 

The spectral responses of the interference filters were obtained 

0 

by calibrating the transmission of each filter from 6000 to 13,000 A. 

This calibration procedure (Whiting et al. 1972) incorporated a Perkin- 
Elmer 210 grating spectrograph and a synchronous amplifier to obtain 
transmission values as low as 10" 5 of peak transmission. Transmission 

O r 

values for wavelengths greater than 13,000 A were smaller than 10 . 

This is primarily due to the very low sensitivity of the photomultiplier 
photocathode surface at these wavelengths. 

A detailed description of these calibrations procedures and radio- 
metric techniques in general is contained in the work by Craig and Davy 
(1969). An additional check suggested by them is to determine the linear 
operating range of each radiometer. All of the radiometers used in the 
experimental work were checked for this effect. This was accomplished 
by selecting a filter whose transmission is well known and placing it in 
and out of the incident beam of light and recording the signals. The 
ratio of these signals gives the transmission of the filter if the radiom- 
eter is operating in the linear regime. As the transmission is repeatedly 
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measured with increasing intensity, an intensity level will be reached 
at which the calculated filter transmission begins to increase. This 
represents the linear limit. All of the experimentally measured signals 
were kept below this limit. 

In order to check the radiometers used in the ultraviolet below 

O 

3000 A for scattered light, two shock-tube tests were performed with 
glass windows rather than quartz windows in the splitter plates. These 
tests demonstrated that scattered light in the radiometers was negligible. 
4.3 Initial Test Conditions 

The test gas which when shock heated produced excited C 2 was 
nominally 85-percent argon and 15-percent acetylene (C 2 H 2 ) . The actual 
gas composition as certified by the supplier was routinely verified by a 
mass spectrometer analysis of a sample of the test gas taken from the 
shock-tube test section. Contaminant concentrations never exceeded SO 
parts/million in any of the analyses. 

The initial pressure, Pj, in the shock tube ranged from 5 to 15 torr. 
However, the majority of the tests were performed at a pressure of 10 
torr. Shock velocities varied from approximately 4.3 to 5.0 km/sec. The 
shock velocity, initial loading conditions, and test gas composition 
were selected in order to maximize the ratio of C 2 particles to other 
species such as C 2 , CH, and particulate carbon. Figure 4.7 is a plot of 
the concentration of some of the carbon-containing species as a function 
of velocity for an initial pressure of 10 torr. These concentrations were 
obtained using the thermochemistry program of Horton and Menard (19b9) 
described in section 3.5. Particulate carbon concentrations for these 
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Velocity, km /sec 


Figure 4.7.- The variation of the concentration of selected carbon- 
containing species as a function of shock velocity. The test gas is 
85 percent Ar and 15 percent C 2 H 2 ; Pi = 10.0 torr. These concen- 
trations were obtained with the thermochemistry program of Horton 
and Menard (1969) . 
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conditions were quite low - less then 10 12 particles/cm 3 and therefore 
do not appear on the plot. As illustrated on the figure, the majority 
of the shock-tube tests at this pressure were made at velocities between 

4.3 and 4.5 km/sec. The C 2 concentration varies slowly over this 
velocity interval. Consequently, a slight error in determining the shock 
velocity will result in only a small error in the C 2 particle 
concentration. 

Note that the species concentrations obtained from the program of 
Horton and Menard (1969) are compatible with those obtained from the 
NASA Lewis program (McBride and Gordon, 1967). 

4.4 Sample Spectrographic and Radiometric Data 
4.4.1 Spectrographic Data 

Photographic survey spectra of the shock-heated gas were obtained 
to verify the purity of the test gas, to demonstrate that C 2 was the 
dominant radiator in the radiometer bandpass, and to identify the C 2 
emission. Sample survey spectra obtained with the spectrographs 
described in section 4.2.3 are displayed in figures 4.8 to 4.10. These 
exposures were each obtained from a single shock-tube test. As men- 
tioned in section 4.2.3, the spectrographs were equipped with explo- 
sively driven capping shutters with a closing time of 35 usee. Therefore 
these spectra represent a time integrated exposure of the shock-heated 
gas. It will be shown below that the shock-tube test time exceeds 
35 usee for all of the test conditions utilized and consequently these 
spectra are from the shock-heated test gas only. 
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Figure 4.8.- Identified survey spectrum of the shock -heated Ar and C 2 H 2 
in the ultraviolet region. This spectrum was obtained with a Huet 
UV 24 spectrograph. The path length was 33.6 cm (double-pass 
viewpath) and the exposure is from a single shock-tube test. The 
spectrum was recorded on a type 103-0 plate. 
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Figure 4.9.- Identified survey spectrum of the shock-heated Ar and C 2 H 2 
in the ultraviolet and visible regions. This spectrum was obtained 
with a Jarrell -Ash model 75-000 spectrograph. The path-length was 
16.8 cm and the exposure is from a single shock-tube test. The 
spectrum was recorded on Tri-X panchromatic film. 
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Figure 4.10.- Identified survey spectrum of the shock-heated Ar and 

C 2 H 2 in the visible and infrared regions. This spectrum was obtained 
with a stigmatic-grating spectrograph which was constructed in- 
house. The path-length was 16.8 cm and the exposure is from a 
single shock-tube test. The spectrum was recorded on high-speed 
infrared film. 
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Six of the band systems of interest - the Freymark, Mul liken, 
Fox-Herzberg, Deslandres-d 'Azambuja, Swan, and Phillips systems - are 
readily identified in these spectra. The Av * 0, -1, -2, and -3 
sequences of the Freymark system appear in figure 4.8. The Mul liken 
system is a headless system (section 3.4.5) and does not have discrete 
structure. Only the Av = 0 sequence appears since it is the strongest 
radiating sequence of the system (see figure 3.9). This sequence in 
figure 4.8 resembles a doublet consisting of diffuse lines, when in 
reality the structure consists of R and P branches. Six sequences - 
Av = +2, +1, 0, -1, -2, -3 - of the Fox-Herzberg system are also iden- 
tifiable in the spectrum in figure 4.8. Additional sequences (Av = -4 
and -5) of the Fox-Herzberg are seen in figure 4.9. The Deslandres- 
d'Azambujt system is represented by four sequences (Av = +2, +1, 0, and 
-1). The Av = +1 and 0 sequences of this system are badly overlapped 
by the CN violet system which appears as a contaminant in the radiating 
test gas. Several sequences of the dominant Swan system appear in the 
spectra in figures 4.9 and 4.10. The Phillips system is identifiable 
by the presence of the Av = +4, +3. and +2 sequences. The only other 
band system of importance to this study is the Ballik-Ramsay system. 

O 

Since the short wavelength limit of this system is about 11,0/0 A, no 
attempt was made to record it photographically. 

In addition to the Cj band systems, several atomic lines of C I 
and Ar I as well as the H a line have been identified. The Messer le- 

O 

Krauss system which radiates in the 3100-3400 A spectral region could 
not be detected in these spectra. 
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These spectra demonstrate that the shock-heated test gas is 
virtually contaminant free. With the exception of a few potassium 
lines (whose origin is probably the shock-tube walls), the only observed 
radiating contaminant is the prolific CN radical. An estimate of the CN 
concentration was performed by integrating the intensity of the CN radia- 
tion at 3883 and 3590 X from a densitometer trace of the spectrographic 
plate and comparing with the absolute intensity measurements of the 
nearby Deslandres-d'Azambuja system. The CN concentration was estimated 
to be approximately 10 13 particles/cm 3 which is four orders of magnitude 
less than the typical C 2 concentration (figure 4.7). 

Spectra similar to those displayed in figures 4.8 to 4.10 were used 
to select the band-pass location and width to assure that C 2 was the 
dominant radiator in the radiometer bandpass . These measurements are 
discussed in the next section. 

4.4.2 Radiometric Data 

Absolute intensities of selected spectral regions of the radiation 
behind the incident shock were measured by calibrated radiometers. 

These radiometers (section 4.2.4) consisted of dispersing elements and 
photomultiplier tube detectors. The output of each detector was ampli- 
fied, displayed on an oscilloscope, and photographed. Sample oscillo- 
grams of the radiation signal for each of the band systems measured in 
this study are presented in this section. 

Oscillograms from a radiometer whose bandpass was centered at 

O 

5165 A (the Lv = 0 sequence of the Swan system) are displayed in fig- 
ure 4.11. The oscillograms in figures 4.11 (a) and (b) were obtained 
from tests in which the initial pressure was 5.0 and 10.0 torr. 
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Figure 4.11.- Sample oscillograms from the radiometer whose bandpass 
was centered at 5165 A (the Av = 0 sequence of the Swan system); 
(a) Pj - 5.0 torr, (h) P] = 10.0 torr. 
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respect ively. The initial sharp rise in the traces represents the 
passage of the luminous shock front past the slit. In figure 4.11(a), 
an initial overshoot attributed to nonequilibrium processes is followed 
by a rapid relaxation to an essentially constant level of intensity. 

This level indicates that the shock-heated gas had equilibrated. 

The nonequilibrium overshoot vanished as the pressure was increased 
to 10 torr (fig. 4.11(b)). The absence of this phenomenon indicated 
that the time required for the gas to equilibrate was shorter than the 
rise time of the instrumentation. The radiometer electronics (photo- 
tube and amplifier) had a rise time of 40 nsec and the signal from this 
unit was displayed on an oscilloscope which had a rise time of 25 nsec. 
However, the optical rise time of the radiometer was about 0.11 usee. 
This rise time was due primarily to the time it took the shock wave to 
move across the slit in the splitter plate. 

Additional oscillograms of the radiation signal from the Freymark, 
Mulliken, Fox-Herzberg, Deslandres-d'Azambuja, Phillips, and Ballik- 
Ramsay systems are displayed in figure 4.12. These oscillograms which 
were obtained from tests with 10 torr initial pressure, are similar to 
that of the Swan system and indicate that for the present study the gas 
had reached equilibrium. 

Oscillograms similar to those just presented can also be used to 
determine the approximate facility flow time. For an initial pressure 
of 5.0 torr, oscillograms with slower sweep rates revealed that the 
intensity behind the shock front remained almost constant for about 
80 usee and then began to gradually decrease. This value is somewhat 
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Figure 4.12.- Oscillograms of the radiation signal from the fa) Bal 1 i k-Ramsav system, fb) Phillips 
system, (c) Pesl andres-d * Azambuja system, fd) Fox-llerzberg system. 
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less than the predicted test time of 130 ysec, based on Mirels' (1963) 
laminar boundary- layer theory. This observed difference between the 
experimental and theoretical test time is probably due primarily to 
mixing at the contact surface (figure 4.1). 
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Chapter 5 

Electronic Transition Moments Measurements and Results 

The object of the experimental program described in this thesis was 
to measure the absolute strength of the band systems of the C 2 molecule 
in the spectral region between 0.2 and 1.2 p. The results of these 
measurements are given in this chapter in terms of the electronic transi- 
tion moments squared values, E|R e /ea5|*. These values were obtained 
using the synthetic spectrum analysis technique described in section 3.6. 

The first section of this chapter deals with the spectral location 
of the radiometers and the dominant Av sequence of each electronic 
transition selected for measurement. The absolute strengths of the 
Ballik-Ramsay, Phillips, Swan, Deslandres-d 'Azambuja, Fox-Herzberg, 
Mulliken, and Freymark systems are presented in section 5.2. An error 
analysis leading to an estimate of the accuracy to which the r|R e /ea 0 P 
values have been determined is presented in section 5.3. 

5. 1 Spectral Measurements 

Survey spectra were used to select the strongest radiating Av 
sequences which were free of contaminant radiation for radiometric 
measurement. Densitometer traces of some of these spectra are shown in 
figures 5.1 to 5.3. The sequences selected for measurement are indi- 
cated in the figures by the short horizontal lines. The length of these 
lines represent the approximate width of the bandpass of the radiometer 
assigned to measure a given spectral feature. Table 5.1 contains the 
wavelengths at which the radiometer bandpasses were centered, the band- 
pass of each, and the dominant Av sequence in the bandpass for each 
band system. 
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Figure 5.1.- Densitometer trace of the ultraviolet survey spectrum. The 
small horizontal bars (m) indicate the spectral location and width of 
the radiometers. 
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Figure 5.2.- Densitometer trace of the visible survey spectrum. The small 
horizontal bars C' — ' ) indicate the spectral location and width of the 
radiometers. 
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Figure 5.3.- Densitometer trace of the infrared survey spectrum. The small 
horizontal bars (i — 0 indicate the spectral location and width of the 
rad iometers . 


- 116 - 


Table S.l Radiometer Settings and Information. 


Band system 

band -pass 
centerline, 
A 

Dominant 
Av sequence 
in bandpass 

Band-pass width, 

X 

Swan 

5165 

0 

19 

d 3 w g - a 3 tt u 




Phillips 

12,140 

0 

110 

A l w u - X x Zg 

8750 

+2 

20 

Des 1 andres -d ' Azambu j a 

4090 

-1 

19 

C 1 n g - A 

3400 

♦2 

28 

Fox-Herzberg 
e 3 ir g - a 3 n u 

2775 

-2 

30 

Mulliken 

D h* - X *e + 
u g 

2315 

0 

30 

Freymark 

E H+ - A 

g U 

2218 

-1 

30 

Ballik- Ramsay 

b V - a 3 ir 
g u 

11,970 

♦2 

92 


The radiometers which are described in Chapter 4 were always used 


in pairs with one centered on the bandhead of interest and the other just 
off the bandhead to measure the background radiation. The bandhead-to- 
background intensity ratios were essentially constant for a given band 
system. These ratios did, however, vary from about 8 to 1 for the Swan 
system to approximately 1.4 to 1 for the Freymark system. 
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5.2 Electronic Transition Moment. Results 

The spectral measurements described above were converted to 
electronic transition moments by using the synthetic spectrum data 
analysis technique previously discussed in Chapter 3. As stated, 
a synthetic spectrum was produced by calculating the intensity (equa- 
tion (3.1)) of each rotational line, distributing this intensity in 
the spectrum by a Voigt profile and then summing. This spectrum was 
then folded with the measured instrument profile (Chapter 4) to predict 
a theoretical signal output. The ratio of measured (corrected for 
background radiation) to theoretical signal output yields the 
experimental value of E|R e /ea 0 | 2 . 

The final results of the measurements performed during this work 
are displayed in table 5.2. The E|R e /ea 0 | 2 values in column four are 
the average of the individual electronic transition moments deduced 
from several experimental tests. A complete listing of the L | R e /ea 0 j 2 
values from each test as well as the pertinent thermochemical quantities 
for each test is contained in appendix E. The values of S | R e /ea 0 | 2 
represent an average for a given Av sequence. This is because the 
bandpass often includes contributions from two or more vibrational 
transitions of a particular band. Note that the 2 | R e /ea 0 | 2 values 
listed in table 5.2 are valid for the specific wavelength (or inter- 
nuclear separation) at which they were determined but are not 
necessarily valid for the entire band system. 

Although two sequences of the Deslandres -d'Azambuja and Phillips 
systems were measured, no attempt was made to measure systematically 
the variation of E |R e /ea 0 | 2 with internuclear separation. 
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Table 5.2 The Experimental Values of the Electronic Transition 

Moments . 


C 2 band system 

Av 

X, 1 

r|R e /ea 0 | 2 , A.U. 

Swan, 

d 3 ir - a 3 ir 

g u 

0 

5165 

3.52 ± 0.50 

Phillips, 

A ! ir - X l E + 
u g 

0 

+2 

12,140 

8750 

0.36 ± 0.08 
0.40 ± 0.10 

Deslandres-d'Azambuja, 

C J Tr - A 

g u 

H 

4090 

3400 

0.93 ± 0.18 
0.91 ± 0.18 

Fox-Herzberg, 
e 3 ir - a 3 tt 

g u 

— 

-2 

2775 

0.40 ± 0.10 

Mulliken, 

D 1 I + - X 1 Z + 
u g 



0 

■ i 

2315 

0.13 ± 0.05 

Ficjymark, 

E *E + - A 1 ir 
g u 

-1 

2218 

2.26 ± 1.13 

Ballik-Ramsay, 

b 3 E” - a 3 tt 
g u 

+2 

11,970 

0.65 ± 0.15 


5.3 Error Analysis 

An estimate of the accuracy to which the ljR e /ea 0 | 2 values have 
been determined has also been included in table 5.2. This uncertainty 
was obtained by summing estimates of individual experimental sources of 
error and therefore represents the maximum error possible due to 































-119- 


experimental uncertainties. Individual sources of error with their 
accompanying estimate of accuracy include: 

Shock velocity, ±0.1 percent 

Initial pressure, ±2.0 percent 

Gas composition, ±0.03 percent 

Band-pass position, ±2 A 

Absolute calibration, typically ±5 percent 

Signal deflection on oscillogram, ±4 percent. 

The uncertainity in the absolute calibration was larger (±10 percent) 
for the systems which lie below 2S00 A. This results from the addi- 
tional uncertainties introduced by extending the calibration curves to 

O 

2200 A. The large uncertainity in the Freymark system results from 
this additional error as well as the small bandhead-to-background signal 
ratio. 

No attempt was made to estimate the uncertainties which would 
arise from errors in the spectroscopic or thcrmochemical constants. 

The information contained in appendix D is sufficient to allow a 
rapid re-evaluation of the data for any significant changes in these 
constants. 
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Chapter 6 

Discussion and Conclusions 

In this chapter, the electronic transition moments presented in the 
preceding chapter are compared with the results found in the literature 
and summarized earlier in section 2.S. A critical comparison with both 
experimental and theoretical results are giver, in section 6.1, and the 
final conclusions are presented in section 6.2. 

6.1 Comparison With Other Results 

A comparison of the final results of this study with both 
experimental and theoretical results found in the literature is pre- 
sented below for each band system. Most of the results cited in the 
literature are defined in terms of various absolute intensity parameters 
and care must be used in interpreting them. The original references 
were consulted in order to determine the manner in which the results were 
defined. These results were then converted to z|R e /ea 0 | 2 values by 
using a consistent set of defining equations. The defining relations 
were taken from Schadee (1967) and Tatum (1967) and are reproduced in 
appendix A. 

6.1.1 The Swan System, d 3 ijg - a 3 n u 

The Swan system has been studied extensively in the past to obtain 
electronic transition moment data. Table 6.1 presents a comparison of 
the present experimental value (3.52 ± 0.50) with others in the litera- 
ture. Only measurements of the Av = 0 sequent ' (5165 X) are included. 
The present results agree excellently with the shock-tube measurements 
of Arnold (1968) and Fairbairn (1966). Both of these measurements were 
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Table 6.1 Measurements of the Electronic Transition Moments of the 
Swan System. Only Measurements of the Av » 0 Sequence (5165 X) 
are Included. 


Experimenter 

Method 

Z|R e /ea 0 | 2 , A.U. 

Present results 

Shock tube 

3.52 ± 0.50 

Arnold (1968) 

Shock tube 

3.56 ± 0.50 

Fairbairn (1966) 

Shock tube 

3.36 ± 1.20 

Harrington et al . (1966) 

Shock tube 

2.95 ± 0.92 

Sviridov et al. (1965) 

Shock tube 

2.24 ± 0.81 

Sviridov et al. (1966) 

Shock tube 

2.85 ± 0.20 

Hicks (1957) 

King furnace 

2.64 

Hagen (1963) 

King furnace 

0.71 ± 0.40 

Jeunehomme and Schwenker (1965) 

1 

Laser 

excitation 

0.60 ± 0.02 

Fink and Welge (1967) 

Phase shift 

2.24 ± 0.70 

Wentink et al. (1967) 

| 

Laser 

excitation 

>1.35 


also made behind the incident shock. Arnold used a test gas consisting 
of CO 2 and Ar, while Fairbairn used a variety of test gases which con- 
sisted of one or more of the following: C 2 H 2 , C 2 N 2 , CO, and Ar. Addi- 

tional shock-tube measurements by Harrington et al. (1966), Sviridov 
et al. (1965), and Sv^idov et al. (1966) obtained in the reflected 
shock region are in fair agreement with the present results. Harrington 
et al. used fluoro-carbon as a test gas, while Sviridov et al. used 
pure CO and mixtures of CO with Ar. 
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There exists an almost uniform discrepancy (albeit small) between 
the incident and reflected shock measurements. This discrepancy appears 
to be more than a coincidence and may result from a common error in cal- 
culating the jump conditions across the shock. If this hypothesis were 
true, one would suspect the reflected shock properties since they are 
more complicated to predict than the single jump properties across an 
incident shock. Fairbairn (1966) has suggested that the uncertain bond 
energies in the fluoro-carbon system might be another source of error 
in the data of Harrington et al . 

The data of Hagan (1963) obtained in a King furnace displays an 
extremely large scatter with z|R e /ea 0 | 2 values ranging from 0.32 to 
4.37. This large scatter is suggestive of nonequilibrium conditions 
which is the major disadvantage of King furnaces. The value (0.71) 
recommended by Hagan is clearly too small. Similar measurements in a 
King furnace by Hicks (1957), however, agree reasonably well with the 
shock-tube results. Note that the value of 3.45 reported by Hicks has 
been adjusted to 2.64 to account for the newer values of the heat of 
sublimation of C 2 (Fairbairn, 1966). 

Jeunehomme and Schwenker (1965) and Wentink et al. (1967) have 
used the laser excitation technique to obtain radiative lifetime measure- 
ments of the upper state of the Swan system. Both measurements were 
made in the hot vapor arising from graphite being heated by a laser beam. 
Wentink et al . discuss two major problems of this technique which 
lead them to conclude that their results usually gives only a lower limit 
to the z|R e /ea 0 | 2 value. These problems are: (1) the temperature and 
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pressure behavior of the blowoff vapor is complex and not well under- 
stood; and (2) the high velocity of the blowoff vapor produces a time 
variation of the density of the blowoff vapor in the field of view of 
the instrumentation optics. Fairbairn in discussing the data of Jeunehommc 
and Schwenker suggested the possibility of formation of molecules in the 
upper state during the decay process as an additional source of error. 

It now seems justified to assume that the results of Jeunehomme and 
Schwenker and Wentink et al . in table 6.1 are incorrect. 

The lifetime measurements of Fink and Welge (1967) using the 
phase-shift method agree w A cn the lower bounds of the shock-tube data 
and are thus in fair agreement with the present results. In this work, 
the lower state molecules were excited by electron bombardment of Collg 
at low C 1 0~ 3 torr) pressures. These results are based on the critical 
assumption that the excited state formation cime is very short compared 
to the radiative lifetime. This assumption as Fink and Welge point out 
is not trivial, and represents a potential source of error in the life- 
time measurement data. 

Arnold (1968) concludes that the correct z|R e /ea 0 | 2 value of the 
Av = 0 sequence of the Swan system is probably bounded by the shock- 
tube measurements listed in table 6.1. The results from the present tests 
support that conclusion and indicate that the correct value lies near 
the upper bound of these measurements. 

The results of various theoretical efforts to predict the absolute 
strength of the Swan system are listed in table 6.2. Clementi (I960) 
lias calculated the strength of several of the band systems of C 2 . Me 
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Table 6.2 Theoretical Predictions of the Absolute Strength of the 


Swan System. 


Reference 

E|R e /ea 0 | 2 , A. U. 

Remarks 

dementi (1960) 

1.14 

z|R e /ea 0 | 2 value modified to 
account for the statistical 
weight of the upper level omitted 
from defining-equation 

Shull (1951) 

4.24 

Same as above. 

Stephenson (1951) 

4.06 


Lyddane, Rogers, 
and Roach (1941) 

3.39 


Spindler (1965) 

5.04 

Based on Jeunehomme's lifetime 
measurements of the d 3 ir g state 

Coulson and Lester 
(1955) 

- -- - - - - - ■ | 

33.92 


Results of this 
work 

3.52 ± 0.50 

Experimental measurement of the 
Av = 0 sequence (5165 A) 


used the dipole-moment operator method and obtained a E|R e /ea 0 | 2 value 
of 1.14 for the Swan system. As noted in table 6.2, dementi omitted the 
statistical weight of the upper state from his defining equation. Con- 
sequently, the £ | R e /ea 0 | 2 value credited to him on table 6.2 has been 
modified to account for this factor. A similar correction was also made 
to the results of Shull (1951). Shull, Stephenson (1951), Lyddane et 
al. (1941), and Coulson and Lester (1955) all used the dipole-moment 
operator method. Spindler's (1965) prediction is based on the lifetime 
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measurement of the upper state by Jeunehomme. The early work of Lyddane 
et al., Stephenson, and Shull agrees surprisingly well with the higher 
values reported from shock-tube tests. These theoretical papers, how- 
ever, also contain predictions of the absolute strengths of other band 
systems of C 2 , CN, and N 2 . Experimental measurements of the strengths 
of these band systems indicates that only the predictions of Stephenson 
are consistently accurate. Since these calculations are quite old and 
based on hydrogen-ion approximated wave functions, it is clear that 
there is still much progress to be made in the field of theoretical 
computations of basic molecular properties. 

A new measurement of the relative spectral intensity distribution 
of the Swan system was recently reported by Danylewych (1972) and 
Danylewych and Nicholls (1974). This work contains a valuable plot 
comparing all of the past measurements of the variation of spectral 
intensity of the Swan system. 

6.1.2 The Fox-Herzberg System, e 3 *g - a 3 tt u 

The value of £|R e /ea 0 | 2 for the Av = +2 sequence (2775 A) of the 
Fox-Herzberg system was found to be 0.40 ± 0.10. Sviridov et al. 

(1966), also using a shock tube but obtaining measurements in the 
reflected shock region, placed an upper limit 0.55 on the S | R e /ea 0 1 2 
value. These measurements were made in mixtures of CO and Ar. Went ink 
et al. (1967) using the laser excitation technique placed a lower limit 
of 0.45 on the Av = 0 sequence (2512 A). (See section 6.1.1 for a 
discussion of this work.) The excellent agreement between these three 
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measurements nay be fortuitous since two of then represent limits. The 
results however are compatible. 

dementi’s (1960) prediction of 57.80 is obviously unreliable. This 
value has been reduced by a factor of six to account for the statistical 
weight of the upper level missing from the equation defining the absolute 
strength in his paper. 

6.1.3 The Deslandres-d'Azambuja System, C - A 1 it VL 

Values of 0.91 ± 0.18 and 0.93 ± 0.18 for the Av « +2 (3400 A) and 
-1 (4090 X) sequences, respectively, were obtained \n the present work. 
The sequence at 3852 X was not measured since the Av = 0 sequence of CN 
violet (which was observed as a slight contaminant) overlaps this wave- 
length region. Wentink et al. (1967) quote a lower bound of 0.18 for 
the Av = 0 sequence. 

Stephenson (1951), Shull (1951), and dementi (1960) have all used 
the dipole moment operator method to predict z|R e /ea 0 | 2 values of 
1.35, 3.97, and 1.12, respectively. Stephenson's prediction is greater 
(46 percent) than the experimental value. However, as discussed above, 
his predictions are consistently 20 to 50 percent higher than the 
experimental values. While dementi's prediction agrees quite well with 
the experimental value, it appears to be a coincidence since his pre- 
dictions of the other C 2 band strengths are not supported by experimental 
measurements nor do they show a consistent difference with the experi- 
mental results. Shull's value of 3.97 is too high. Once again the 
results of dementi and Shull have been modified (section 6.1.1) to 
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account for the incorrect definition of the absolute strength used by 
them. 

As stated earlier, a detailed study to determine the variation of 
E |!< e /ea 0 J 2 with internuclear separation, r, was not undertaken. How- 
ever the results of the present tests indicate that the E | R e /ea 0 | 2 
value does not vary with r over the range of r (and A) appropriate 
to the two sequences of the Deslandres-d'Azambuja system that were 
measured. 

6.1.4 The Mulliken System, D - X 1 E + 

^ 8 

A value of 0.13 ± 0.05 was measured for the E | R e /ea 0 | 2 value of 
the Av = 0 sequence (2315 K) of the Mulliken system. Smith (1969) 
measured the radiative lifetime of the same sequence using the phase 
shift method. He produced C 2 spectra by focusing high energy electrons 
on ethylene or benzene. Smith obtained a E | R e /ea 0 [ 2 value of 
0.42 ± 0.05. These lifetime measurements are based on the nontrivial 
assumption that the excited state formation time is very short compared 
to the radiative lifetime. This assumption was not verified by Smith 
and represents a potential source of substantial error. It is interest- 
ing to note also that Smith was unable to detect the Freymark, Fox- 
Herzberg, and Deslandres-d'Azambuja systems contained in the spectral 
region he examined, dementi's (1960) prediction for this system is 
0.78. 

6.1.5 The Phillips System, A 1 ir u - X 

The r|R e /ea 0 |2 value for the Av * 0 (12,140 A) and +2 (8750 A) 
sequences was determined in the present tests to be 0.36 ± 0.08 and 
0.40 ± 0.10, respectively. These measurements were obtained using the 
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two types of radiometers described in section 4.2.4. The Av = *2 
sequence was measured with the 1/4-m monochromator type. Hicks (1957) 
has also measured the Av = +2 sequence using a King furnace as a 
source. His measurement was not an absolute determination per se but 
was taken relative to the (1,0) state of the Swan system. This measure- 
ment yields a E|R e /ea 0 | 2 value of 0.68 with a factor of three uncer- 
tainty. Considering the gross uncertainty in Hicks' result, it would 
be premature to form any conclusion pertaining to the absolute strength 
of this sequence. However, the two sets of measurements do agree 
satisfactorily, dementi's (1960) prediction of the E|R e /ea Q | 2 value 
for this system is 0.13. 

Based on the measurements of the Av = 0 and +2 sequences, the 
zjRg/eaQl 2 of this system appears to have a weak dependence on the 
internuclear separation. 

6.1.6 The Ballik-Ramsay System, b 3 E~ - a 3 ir u 

A E|R e /ea 0 | 2 value of 0.65 ± 0.15 was measured for the Av = +2 
sequence (11970 X) of this system. There are no other existing measure- 
ments of this system known to the author. Clementi (1960) has predicted 
a value of 1.21. 

6.1.7 The Freymark System, E *E* - A *ir u 

The E | R e /ea 0 | 2 value for the Av = -1 sequence (2218 A) was 
determined to be 2.26 ± 1.13. There are no other measurements or 
predictions of the strength of this system known to the author. 

A summary of the results presented in sections 6.1.2 to 6,1.7 is 
given in table 6.3. The results of the present effort (excluding the 



Table 6.3 Experimental and Theoretical Results of the Absolute Strength of the Band Systems (Exclud- 


ing the Swan System) of the C 2 Molecule. 


Band system 

Results in literature, 
E|R e /ea 0 | 2 , A.U. 

Results of 
this work 
E|R e /ea 0 | 2 , A.U., 
[Av/A(£)J 

Theoretical 

Experimental 

[av/a(A)J 

Bal lik-Ramsay, 
b 3 Zg - a 3 ir u 

1 .21 -Clement i (1960) 


0.65 ± 0.15 
(+2/11,970) 

Phillips, 

A - X x Eg 

0.13-Clementi (1960) 

0.68-Hicks (1957) 
(+2/8750) 

0.40 ± 0.10 
(+2/87S0) 

Des landres -d 9 Azambu j a , 
C l w g - A 

1.12-Clementi (1960) 
3.97-Shull (1951) 

1 .35-Stephenson (1951) 

>0. 18-Wentink et al. (1967) 
(0/38S2) 

0.91 ± 0.18 
(+2/3400) 

Fox-Herzberg, 
e 3 * g - a 3 » u 

57.80-Clementi (1960) 

<0.55-Sviridov et al. (1965) 
>0.45-Wentink et al. (1967) 
(0/2512) 

0.40 ± 0.10 
(-2/2775) 

Mul liken, 

D K - X '*1 

0.78-Clementi (1960) 

0.42-Smith (1969) 
(0/2315) 

0.13 ± 0.05 
(0/2315) 

Freymark, 

E 1 Eg - A *ir u 



2.26 + 1.13 
(-1/2218) 
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Swan system) are compared with the experimental and theoretical results 
found in the literature. 

6.2 Conclusions 

A systematic study of the absolute strength of seven of the band 
systems of the Cj molecule in the spectral region 0.2 to 1.2 u has been 
made. Experimental measurements of the strength of the Ballik-Rarasay, 
Phillips, Swan, Deslandres-d’Azambuja, Fox-Herzberg, Mulliken, and 
Freymark systems were performed. 

A combustion-driven shock tube was used to heat the test gas of 
85-percent Ar and 15-percent C 2 H 2 to temperatures of approximately 6000 6 K. 
Measurements of the absolute intensity of the equilibrium radiation in 
selected spectral regions were obtained with spectrally and absolutely 
calibrated narrow pass-band radiometers. The measurements were made 
behind the incident shock. 

The absolute intensity measurements were converted to electronic 
transition moments by a synthetic spectrum analysis. The values deduced 
from this work are as follows: 


System 

A (A) 

Ballik-Ramsay (Av * +2) 

11,970 

Phillips (Av » 0) 

12,140 

(Av * +2) 

8750 

Swan (Av * 0) 

5165 

Deslandres-d'Azambuja (Av * +2) 

3400 

(Av * -1) 

4090 


z|R e /ea 0 | 2 , A.U. 
0.65 ± 0.15 
0.36 ± 0.08 
0.40 ± 0.10 
3.52 ± 0.50 
0.91 ± 0.18 
0.93 ± 0.18 
0.40 ± 0.10 


Fox-Herzberg (Av * +2) 


2775 



System 

Mul liken (Av - 0) 
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(X) 


(A) E|R e /ea 0 | 2 , A.U. 


2315 0.13 ± 0.0S 

Freymark (Av » -1) 2218 2.26 ± 1.13 

These results are valid for the specific sequence (or wavelength) at 
which they were determined but are not necessarily valid for the entire 
band system. 
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APPENDIX A 

Absolute Intensity Relationships 
The relationships between the various absolute intensity 
parameters are given herein. Complete descriptions of these relations 
can be found in the papers of Schadee (1967) and Tatum (1967). 

A wavelength-dependent electronic absorption f number, f e ^(A), 
is related to the E|R e /ea 0 | 2 value by 

It c 

e 


8 n 2 m<i c E 1 R e / ea 0 [ 2 


f AX) = 

e * 3he 2 A (2 - « 0 >A „)(2S» ♦ 1) 


, dimensionless. 


The band oscillator strength, f yl „ (written for absorption), and 
the band Einstein coefficient, A v , v „, are related to the electronic 

transition moment by 

8 n 2 m c q , E |R (r . „) 1 2 

e M v’v" 1 e v v'v" 1 , . . 

f , „ = — y , dimensionless 

V V 3he 2 v'v" (2 - S q A „) (2S" + 1 ) 


and 


A _ 64 TT* t V 3 

v'v" 7 , 3 ^v'v" 


E l R e^Vv» ) 


- f sec 


-1 


(2 - « oA ,)(2S’ . 1) 


The band strength, S t lf , is related to the E|R e /ea 0 j 2 through 


v'v 


^ v iyf» ~ ^ Re/ ea o 1 A.U. 


The mean radiative lifetime t , (written for the upper state) is 


given by 


64 tt 


T , £ v 3 q . E|R (r , „)| 2 , sec -1 

T v' 3hc 3 (2 - 6 q a ,)(2S* + 1) ~ 


yi A v »v»» sec 1 
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APPENDIX B 
Band Origins 

The band origins of the seven C 2 bands systems of interest in this 
work are tabulated in this appendix (tables B-l through B-7) . In 
section 3.1.2, the band origin, v , was defined as: 

v (v , ,v")(cm- 1 ) = (T ' - T ") + [G'(V) - G"(v")] = 10 8 /X (A) 
w C “ o 

where T e and G(v) are the electronic and vibrational term energies, 
respectively. The electronic term energies are listed n table 2.1 and 
the vibrational term energies were calculated from equation (3.11) using 
the vibrational constants also listed in table 2.1. These band origins 

O 

are calculated in vacuo and are listed in A units. 


TABLE B-l. BAND ORIGINS (A, IN VACUO) OF THE FREYMARK SYSTEM 



TABLE B-2. BAND ORIGINS (A, IN VACUO) OF THE MULLIKEN SYSTEM 
































































TABLE B-3. BAND ORIGINS (l, IN VACUO) OF THE FOX-HERZBERG SYSTEM 



Hi 

0 


2 

3 

4 

i 

5 

6 

7 

8 

9 

10 

0 

2512.15 

2618.59 

2732.70 

2855.31 

FEEMffrl 

3129.91 

3284.23 





1 

2448.36 

2549.35 

2657.39 



3031.50 


wEsam 




2 

2390.77 

2486.97 

2589.68 



2943.70 

3079.81 

3226.68 

3385.59 



3 

2337.82 

2429.72 

2527.66 

2632.21 

2744.01 

2863.83 

2992.49 

3130.97 

3280.37 

3441.97 


4 

2288.19 

2376.16 



2675.89 

2789.70 

2911.65 

3042.58 

3183.48 

3335.46 

3499.81 

5 


2325.05 

2414.58 

2509.80 

2611.25 

2719.52 

2835.28 

2959.29 

3092.41 

3235.62 

3390.06 

6 

| 











7 
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TABLE B-4. BAND ORIGINS (A, IN VACUO) OF THE DESLANDRES-d ' AZAMBUJA SYSTEM 



0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0 

3850.62“ 

4100.77 

4381.01 

4697.00 

5055.89 

5466.84 






1 

3605.65 

3824.08 

4066.66 

4337.53 

4641.81 . 

4985.91 

5378.03 





2 

3397.79 

3591.09 

3804.19 

4040.21 

4302.94 

4597.04 

4928.34 

5304.19 




3 

3222.11 

3395.42 

3585.32 

3794.22 

4025.01 

4281 .22 

4567.15 

4888.13 

5250.85 



mmm 

3074.63 


Knsnoi 


iiMrltll 

IHsflBB 

4276.40 

4556.56 

4870.17 

5223.39 



2952.23 ! 

3097.07 

3254.29 

3425.48 

3612.48 

3817.53 

4043.27 

4292.79 

4570.03 

4879.68 

5227.56 

6 

... _ j 

2987.50 

3133.53 

3291.94 

3464.28 

3652.41 

3858.49 

4085.12 

4335.40 

4613.10 


7 





3348 . 76 

3524.23 

3715.72 

3925.43 

EUdA&I 

4410.49 

4692.76 

8 





E3SSS 

3429.56 

3610.63 

E£TS!3E!9 

4024.96 

WESSk 9 

4526.38 

9 





3205.97 

3366.44 

3540.75 

3730.67 

3938.30 

4166.12 

4417.09 

10 








3505.27 

3691.30 

3894.45 

4117.08 











































































































































TABLE B-5. BAND ORIGINS (A, IN VACUO) OF THE SWAN SYSTEM 



v \v" 

v * N 

0 

1 

1 

2 

3 

4 

5 

6 i 

l 7 

8 

9 

1 

10 

0 

5160.37 

EgaBKH 

■aoEaiai 


7664.92 

8679.21 

■ ■ 

B B 

BiWWf 

B B 

BBbHBi 

1 

4732.13 

5124.49 

5580.53 

wsnmn 

6756.69 

7532.68 

8493.16 





2 

4376.68 

4710.23 

5092 . 76 

5535.76 


6670.34 

7412.66 

8324.69 




3 

4077.61 

4365.64 

4692.31 

5065.82 

5496.87 

5999.70 

6593.61 

7305.55 




4 

3823. IS 


4358.49 

4678.93 

5044.28 

5464.55 

5952.93 

6527.21 

7211.96 

8042.10 


5 

3604.63 

3827.88 

4076.74 

4355.76 

4670.69 


5439.47 

5914.99 

6471.84 

7132.54 

7928.76 

6 


3615.31 

3836.49 

4082.60 

4358.02 

4668.20 

5020.02 

5422.33 

5886.64 

6428.25 

7067.94 . 

warn 




3849.47 

4093.39 

4365.86 

4672.09 

5018.64 


5868.62 

6397.20' 

ES 




3648.92 

3867.36 


4379.92 

4683.08 








i 

3672.90 

3890.79 

4132.15 

4400.92 

4701.92 

5041.19 

5426.33 

10 




1 


3702.46 

3920.37 

4161.50 

4429.64 

4729.50 

5066.90 


TABLE B-6. BAND ORIGINS (A, IN VACUO) OF THE PHILLIPS SYSTEM 


\v" 

0 

1 

— 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0 

112094.60 

15526.30 

21545.02 









1 

10149.90 

12461.30 

16062.70 

22443.55 








2 

8762.53 

10433.23 

12844.36 

1 6623.61 

23387.17 







3 , 

7723.30 

8992.51 

10728.31 

13243.02 

17207.39 







4 

6916.07 

7916.65 

9231.59 

11034.62 

13656.17 

17811 .56 






5 

6271.20 

7082.94 

8117.42 

9479.39 

11351.36 

14082.23 






6 

5744.41 

6418.17 

7256.09 

8325.33 

9735.35 

11677.49 






7 

5306.16 

5875.94 

6570.60 

7435.32 

8539.98 

9998.73 

12011.65 





8 

4936.03 

5425.43 

6012.33 

6728.34 

7620.32 

8760.83 

10268.59 

12352.06 




9 

4619.43 

5045.34 

5549.08 

6153.46 

6891.17 

7810.70 

8987.19 

10543.71 

12696.61 

J 


02 

4345.65 

4720.54 

5158.68 

5677.04 

6299.17 

7058 . 78 

8005.93 

9218.19 

10822.61 

13042.64 
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TABLE B-7. BAND ORIGINS fA. IN VACUO) OF THE BALLIK-RAMSAY SYSTEM 


SI 

0 

1 

2 

3 

4 

5 

6 

7 

, 

8 

9 

i 

10 


IBBBBa 

EEMBBS 



■ ■ 




■ 

■ 

MjlPIBHI 

i 












2 

117SS.34 












3 

10090.18 

12058.92 




■ 

■ ■ 

■ 


Bpili Hi 


4 

8855.42 

10336.44 

12376.48 









S 

7903.53 

9062.43 

10593.34 

12708.80 








6 1 

7147.50 

WmBSn 

9277.96 

10861 .42 

13056.68 







7 




9502.42 

11141.26 







8 





9736.24 

11433.47 






1_2 






9979.89 

11738.69 





10 





i . 


10233.85 

12057.60 




li 





! """ — 



10498.62 

12390.92 
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APPENDIX C 

Franck-Condon Factors 

The Franck-Condon factors used in the data reduction in this study 
are listed in this appendix. These factors were defined in section 3.2 
by 

Vv" " dr l 2 

where the iJj's represent the vibrational wave functions. The factors 
listed herein are based on Klein-Dunham potentials and were taken from 
the work of McCallum, Jarmain, and Nicholls (1970). 


TABLE C-l. FRANCK-CONDON FACTORS FOR THE BALLIK-RAMSAY SYSTEM 




0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 


6. 3400E-01 

3. 0600E-01 

5.4900H-02 1 

4 . 7600E-03 

2. 0500E-04 

2.9400E-06 


1 . 4000E-08 

2.9400E-11 

3.4900E-09 

9.3000E-10 

1 

2.6900E-01 

1 .8700E-01 

3. 9500E-01 

1 .3100E-01 

1 . 6400E-02 

9. 1800E-04 


5.4100E-08 

8 . 3600E-09 

2 . 8100E-10 

8 .4300E-10 

2 

7.4300E-02 

2 . 9500E-01 

2. 2000E-02 

3. 6400E-01 

2. 0600E-01 

3. 5800E-02 


6.4400E-05 

2 .4100E-07 

2.0100E-08 


3 

1 . 7400E-02 

1.4S00E-01 

2 , 1900E-01 

4 . 0900E-03 

2.8000E-01 

2.6700E-01 


5. 2700E-03 

1 . 6400E-04 

9. 8500E-07 


4 

3 . 7900E-03 

4 8800E-02 

1 .8100E-01 

1 .2500E-01 

4.6600E-02 

1 .8500E-01 


9.3400E-02 

9 .6200E-03 

3.4200E-04 

QtTTj ? 

5 

7. 9800E-04 

1 . 3800E-02 

8.4000E-02 

1 . 8000E-01 

S. 0900E-02 



3 . 2300E-01 

1 . 2800E-01 

1 . 5700E-02 

BBS 3SS 

6 

1 . 6200E-04 

3 . 5500E-03 

3 . 0000E-02 

1 . 1 300E-01 

1 .5100E-01 



4 . 3800E-02 

3.2000E-01 

1.6400E-01 

2.3700E-02 

7 

3 . 1900E-05 

8 . 5500E-04 

9 . 2000E-03 

5 . 0000E-02 




1 . 5700E-01 

1 .0100E-02 

2.9700E-01 


8 

6. 7400E-06 

2. 0600E-04 

2 . 6700E-03 

1 . 8S00E-02 




1.2500E-02 

1 . 5200E-01 

1 . 3800E-07 

2.6300E-01 

9 

1 . 5000E-06 | 

4 . 9400E-05 

7 . 3500E-04 

6 . 1S00E-03 

3.0600E-02 



3. 0700E-02 

3. 6700E-02 

1 . 3100E-01 

nreiiiaip 

10 

3 . 5200E-07 

; 1 . 2100E-05 

2.0000E-04 

1 . 9300E-03 

1 . 1700E-02 



9.4700E-02 

8.4000E-03 

6 . 3100E-02 

DEzEiSS 


TABLE C-2. FRANCK-CONDON FACTORS FOR THE PHILLIPS SYSTEM 



0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0 

4 . 1200E-01 

3. 9800E-01 

1 . 5500E-01 

3. 1 100E-02 

3.4300E-03 

2.0500E-04 

5. 8900E-06 

6. 0400E-08 

2.8200E-11 

6. 9300E-12 

' sr sk 

i 

3. 3100E-01 

5. S300E-03 

2.8900E-01 

2 . 71 00E-01 

8.8600E-02 

1 . 3600E-02 

1 . 0300E-03 

3.5600E-0S 

4.0200E-07 




1.6200E-01 

1 . 7200E-01 

S. 8100E-02 

1 . 1 500E-01 

3. 0200E-01 

1 . 5600E-01 

3. 2200E-02 

3.0400E-03 

1.2200E-04 




6.3100E-02 

1 .9900E-01 

3. ] 700E-02 

1.4300E-01 

1.6700E-02 

2.6300E-01 

2 . 1600E-01 

5 . 9100E-02 

6.7900E-03 


tjK ! W 


2. 1600E-02 

1 . 2600E-01 

1 . 3800E-01 

1 . 2300E-03 

1.5600E-01 

2. 7900E-03 

1 . 8900E-01 

2 . 5900E-01 

9.2600E-02 



S 

6.8800E-03 

6. 0200E-02 

1 .4500E-01 

5. 7700E-02 

3.6800E-02 

1.1400E-01 

3. 7300E-02 

1 . 1 100E-01 

2.7800E-01 


2.1300E-02 

6 

2.0900E-03 

2.4600E-02 

9.6100E-02 

1 . 2000E-01 

8 .8S00E-03 

8 . 1500E-02 

5 .7300E-02 

8.3000E-02 

4 .8700E-02 


1 .6800E-01 

7 

6. 1S00E-04 

9.1700E-03 

4 . 9900E-02 

1 . 1 200E-01 

7.3100E-02 

1.6900E-03 

1 . 0300E-01 

1.5S00E-02 

1 . 1600E-01 


2. 5300E-01 

8 

1 . 7800E-04 

3. 2100E-03 

2 . 2500E-02 

7.4200E-02 

1.0400E-01 

2. 9200E-02 

2.2600E-02 

9. 6100E-02 

1 . 1100E-04 i 



9 

5.1400E-0SI 

1 1 .0800E-03 

9.2400E-03 

4.0200E-02 

8.8700E-02 

7 .7700E-02 

4 . 1200E-03 

5.1200E-02 

6.9600E-02 ' 


1.1400E-01 

10 

1 .4900E-05 ! 

3. 5600E-04 

3.5800E-03 

1 .9300E-02 

5. 7900E-02 

8.8500E-02 

4 .4600E-02 

1 . 3700E-03 

7 . IS00E-02 




-6£I“ 







































to vo 


TABLE C-3. FRANCK-CONDON FACTORS FOR THE SWAN SYSTEM 




SSI 

0 

1 


7“" 

L _ 


— 

5 

6 

7 

8 

9 

10 

0 

7. 2400E-01 





2 . 7900E-04 

5.5700E-05 

1 . 2400E-05 

3.1500E-06 

9. 1400E-07 

2.9900E-0? 

1 

2 ,460Gc-01 

3. 3000E-01 




6. 2500E-03 

1 . 4400E-03 

3.4500E-04 

8 .9900E-05 

2.5900E-05 

8 . 3900E-06 







5 . 1 700E-02 

1.S300E-02 

4 . 2700E-03 

1 . 2100E-03 

3.6200E-04 

1 . 1800E-04 


1 . 3500E-03 

7.4700E-02 


2.230^5-02 


1 .6700E-01 

7 . 7600E-02 

2.8400E-02 

9.4600E-03 

3.11 OOE-03 

1 . 0600E-03 


1 . 2900E-05 





1 .2900F-01 

1.6500E-0I 

9. 8800E-02 

4.4100E-02 

1 . 7300E-02 

6.5500E-03 

5 

3.1400E-07 

3.9900E-05 


KillS 


1 . 4800E-02 

6 . 7200E-02 

1 .4200E-0J 

1 . 1000E-01 

5 . 9500E-02 

2 . 7300E-02 

6 

5. 7100E-08 

3 .2600E-06 




3 . 1 500E-01 

l . 2500E-02 

2.4700E-02 

1 . 0900E-01 

1 . 1000E-01 

7 . 1600E-02 

7 

1.8800E-11 





2. 8400E- 01 

2 . 6500E-01 

6 . S000E-02 

3 . 5800E-03 

7 . 1800E-02 

9.8400E-02 

8 


1 .4300E-09 




3 . 1 800E-02 

3 . 2800E-01 

2 . 2200E-01 

8.3400E-02 

7.5400E-04 

3.9500E-02 

9 

2.3300E-11 





7.2800E-05 

3 . 6800E-02 

3 . 6900E-01 

1 . 9100E-01 

I8.6100E-02 

1.0300E-02 

10 


1 .4200E-1 1 


mmsamwm 

4 . 6500E-04 

| 7.0S00E-04 1 

3.7800E-02 

4.0500E-01 

1 1 . 7300E-01 

7 . 7600E-02 


TABLE C-4 . FRANCK-CONDON FACTORS FOR THE DESLANDRES-d 'AZAMBUJA SYSTEM 


6S 

0 

1 

2 

3 


_ 5 

‘ . 

— 

7 

8 

9 

10 

0 



njrepT»T®w| 





7 . 4100E-06 

2.8100E-06 

1.4000E-06 

7.6100E-07 

1 








2 . 2900E-04 

7. 1800E-05 

2.9600E-05 

1 . S200E-05 

2 

2. 7800E-02 

171 


a 



IiB!n Sp 


1 . 0000E-03 

3. 7900E-04 

1 . 7500E-04 

3 

1 . 08 00E-03 

6.8000E-02 


4 . 8200E-02 



Hi!i» ?! wfl 


8.4200E-03 

3. 1600E-03 

1 . 3800E-03 

4 

2. S400E-06 

2.49005-03 





1 .5800E-01 

8.8100E-02 

4 . 0800E-02 

1 . 7100E-02 

7 . 8300E- 03 

5 

1.0500E-05 


HHMJgMHBMiiaro 

4 . 61 00E-01 


gngf;f» iHni 1 

1 . 3300E-01 (8.9100E-02 

5.5400E-02 

2 . 7900E-G2 

6 

1 .8200E-07 







3 . 7100E-03 11 . 0100E-01 

6.6900E-02 

6. 3000E-02 

7 

3.0600E-07 







7 . S800E-03 

jl . 3800E-02 

7 .8100E-02 

2 . S900E-02 

8 

4.4700E-09 







EEZE3SZ1 

12.0700E-02 

9.1300E-02 

7.4800E-02 

9 

3.5300E-08 



1 

8 . 5S00E-04 


2 . 1 700E-04 

jl . 7600E-01 

1 . 5100E-02 

1 . 5700E-01 


6.0100E-09 






5. 9400E-02 !4 .5300E-02 

;1 . 8100E-03 

|l . 0500E-02 


-OH" 























TABLE C-5 


FRANCK-CONDON FACTORS FOR THE FOX-HERZBERG SYSTEM 


Si 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0 

wwtmar«i 

1 . 7100E-02 

6. 1200E-62 

1.3300E-01 

1.9800E-01 

2 . 1400E-01 

1 . 7500E-01 

1.U00E-01 

S.5800E-02 

2 . 2700E-02 

7.6100E- 03 

1 

8.9200E-03 

5 . 0900E-02 

1 . 1700E-01 

1 . 3300E-01 

6. 5300E-02 

2.0900E-03 

3 . 3400E-02 

1 . 2200E-01 

1 . 7000E-01 

1.4S00E-01 

8 . 9000E-02 


2.0200E-02 

8 .2600E-02 

1.1300E-01 

4 . 5900E-02 

8.9100E-04 

6 . 1300E-02 

9 .4300E-02 

3 . 3300E-02 

2.1600E-03 

6.8300E-02 

1 .4400E-01 


3. 3500E-02 

9 . 6400E-02 

6.8300E-02 

7. 3600E-04 

4 . 6500E-02 

7 . 1 200E-02 

8 . 4300E-03 

2 . 5400E-02 

8.6400E-02 

5 . 1 100E-02 

1 .4600E-04 


4 .5300E-02 

8.9400E-02 

2.3500E-02 

1.4900E-02 

6.6500E-02 

1.4300E-02 

2 . 0100E-02 

6 .8100E-02 

1 . 7200E-02 

1 .4500E-02 

7.8800E-02 

5 

4.8800E-02 

6.4100E-02 

2.1400E-03 

3. 7600E-02 

3 . 6200E-02 

2 . 1000E-03 

5 . 0300E-02 

1 . 9100E-02 

1 . 1200E-02 

5.8300E-02 

1 . 8400E-02 

6 

4.4000E-02 

6.8400E-02 

4 . 9700E-03 

3 . 1 000E-02 

3. 3300E-02 

3. 2800E-03 

4 . 9000E-02 

7 . 5300E-03 

2.9400E-02 

4 . 7300E-02 

7.7100E-05 

7 

5. 7200E-02 

S.8300E-02 

3.9000E-05 

4 . 4400E-02 

9 . 9400E-03 

2 .4900E-02 

3. 1 100E-02 

5 . 3800E-03 

4 . 9200E-02 

3.5300E-03 

4. 0300E-02 

8 

6.4700E-02 

4 . 1 200E-02 

6 . 0400E-03 

4 . 2200E-02 

l . 7000E-06 

4 . 0200E-02 

4 . 5800E-03 

3. 2800E-02 

2 .0200E-02 

1 . 7000E-02 

4 . 3800E-02 

9 

6.6600E-02 

2 . 2900E-02 

1 . 7800E-02 

2.8100E-02 

7 . 3400E-03 

3 . 3600E-02 

2 . 1 700E-03 

3.9900E-02 

2 .4900E-06 

4.4200E-02 

4 .6900E-03 

10 

6. 7400E-02 

9.5700E-02 

2.9700E-02 

1 . 2600E-02 

2 . 1 700E-02 

1 . S700E-02 

1 . 8200E-02 

2.07G0E-02 

1.4800E-02 

2.9100E-02 

9.4600E-03 


TABLE C-6. FRANCK-CONDON FACTORS FOR THE MULL I KEN SYSTEM 


si 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 


mbm«h 


1.6200E-04 

3.6400E-06 

1 .7600E-08 

8.3600E-09 

5 .4300E-09 

1 .5900E-09 

3.3200E-10 

5.9100E-H 

1.0100E-11 

1 

2.S900E-03 

9.9300E-01 

3.9400E-03 

4 . 2100E-04 

1 . 1800E-05 

4.9000E-08 

4.9700E-08 

3. 2300E-08 

1 . 0600E-08 

2.7900E-09 

5.5200E-10 

2 

9.1400E-05 

4. 1800E-03 

9.9000E-01 

4 . 5100E-03 

7.2800E-04 

2. 4000E-05 

8.3400E-08 

1 . 7300E-07 

1 . 1900E-07 

4 . 0100E-08 

1 . 0300E-08 

3 

3.7500E-08 

2.6500E-04 

4.9200E-03 

9.8900E-01 

4.4300E-03 

1.0500E-03 

3.8700E-05 

9.3100E-08 

4 . 3900E-07 

3.1300E-07 

1 . 1100E-07 

4 

1.0700E-07 

3. 1200E-07 

5.0900E-04 

4.9600E-03 

9 . 8900E-01 

3.9000E-03 

1 . 3600E-03 

5 .4600E-05 

7 . 2300E-08 

9. 1500E-07 

6.6100E-07 

5 

1.9400E-08 

4.6600E-07 

1 . 2700E-06 

8 .0400E-04 

4.4700E-03 

9.9000E-01 

3.0900E-03 

1 .6600E-03 

7.0400E-05 

4 .0600E-08 

1.6200E-06 


2. 2500E-09 

1 . 1600E-07 

1 . 2200E-06 

3.6800E-06 

1 . 1300E-03 

3. 6200E-03 

9. 9100E-01 

2.1600E-03 

1 . 9200E-03 

8 .6200E-05 

2.1500E-08 

mm 

2.8800E-10 

1.6700E-08 

3.9200E-07 

2.4200E-06 

8.6200E-06 

1 .4700E-03 

2 .S900E-03 

9.9200E-01 

1.2700E-03 

2.1600E-03 

1 .0300E-04 

■s 

3.6500E-11 

2.2200E-09 

6 . S000E-08 

9.9100E-07 

4 . 0200E-06 

1 . 7600E-05 

1.8000E-03 

1 . 5700E-03 

9.9400E-01 

5.S700E-04 

2.3800E-03 

9 

3.9400E-12 

2.8100E-10 

9.6200E-09 

1.9300E-07 

2. 07 OUE- 06 

5. 7? 00E- 06 

3. 240OE-0S 

2. 1000E-03 

7 , 1700E-04 

9.9400E-01 

1 . 2000E-04 

10 

6.2200E-13 

3.4600E-11 

1 . 5100E-09 

3. 1 100E-08 

4 . 5800E-07 

3. 7500E-06I 

7 . 2300E-06 

5 .6200E-05 

2.3700E-03 

1.7500E-04 

9.9400E-01 
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S.8100E 

3.7100E 

4.6900E 

1.5200E 

2.9400E 

1.5600B 

2.8100E 


E-0 

01 1.6200E-01 

02 4.8200E-02 

05 8.4000E-02 
04 S.0800E-04 

06 1.0700E-03 
06 6.9300E-05 


2 

3 

1 .0200E-01 
2.2200E-01 
4.6800E-02 
S.2600E-01 
9.5600E-02 
S.4900E-03 
1.6300E-03 

3.2806E-02 
1.4400E-01 
1.4000E-01 
1 .9500E-02 
5 .6300E-01 
7.7400E-02 
2.0800E-02 


9.9600B- 


03 2.6700E-03 
02 2.4300E-02 

01 8.1700E-02 

02 1.2600E-01 
02 4.0500E-02 

01 3.3900E-02 

02 6.0800E-01 
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APPENDIX D 

Thermochemical Calculation Data 

This appendix contains a list of the species and significant constants 
required as input data in the thermochemical calculations of mixtures of 
argon and acetylene. The constants for atomic species considered in the 
thermochemical calculations are listed in table D-l, while those for the 
diatomic and polyatomic molecules are given on tables D-2 and D-3, 
respectively. 
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TABLE 0-1 . ATOMIC SPECIES DATA USED IN THERMOCHEMICAL CALCULATIONS 


Species 

Number of 
energy levels 
included 

Energy level 
reference 

Heat of 
formation, 
hi 0 , kcal/mole 

C 

41 

Moore 

(1949) 

169.99 

C + 

36 



429. S4 

C++ 

20 



991.69 

C* 

SI 



144.00 

H 

6 



51.62 

H + 

0 



365.03 

H" 

16 

Moore 

(1949) 

344.00 

Ar 

S9 



0.0 

kr* 

28 



363.34 

At** 

41 



1000.28 

e" 

0 



0.0 


TABLE D-2. DIATOMIC MOLECULAR DATA USED IN TOERMOCHEMICAL CALCULATIONS 


Species 

Energy levels 
included 

Energy level 
reference 

Heat of 
formation, 
hi 0 , kcal/mole 

C 2 

X j eJ, a 3 n u , b H' g , 
A c 3 E*, d 3 ir g , 

C l ir g , e 3 Tg, D 1 l*, 

E % 

Ballik and 
Ramsay (1963) 

196.00 

h 2 

X l Z+ 

g 

Herzberg 

(1950) 

0.0 

CH 

X 2 ir, A 2 A, 

B 2 E‘, C 2 E + 

Herzberg 

(1950) 

141.60 
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TABLE D-3. POLYATOMIC MOLECULAR DATA USED IN THERMOCHEMICAL CALCULATIONS ] 

i 

; 


Species 

Ground 

state 

Number of 
vibrational 
frequencies 
included 

Reference 

Heat of 
formation, 
h^°, kcal/mole 

C 3 

h* 

g 

4 

McBride et al. (1963) 

188.00 

c 2 h 

2* 

4 

JANAF Tables (1971) 

116.70 

ch 2 

3 Z" 

g 

4 

McBride et al. (1963) 

95.00 

ch 3 

2 A2 

6 

McBride et al. (1963) 

33.40 

c 2 h 2 

g 

7 

Herzberg (1966) 

54.33 
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APPENDIX E 

Tabulation of Individual Test Conditions and Data 
This appendix contains a complete list of the pertinent thermo- 
chemical quantities, the measured absolute value of the spectral emission 
in the radiometer bandpass, and the l|R e /ea 0 | 2 value deduced for each 
test. The average of the absolute strength factors, r|R e /ea 0 | 2 , given in 
tables E-l to E-7 for the Swan, Fox-Herzberg, Deslandres-d'Azambuja, 
Mulliken, Phillips, Freymark, and Ballik-Ramsay systems, respectively, 
were reported earlier in table 5.2. 

The critical radiometer settings (band-pass center location and 
width) for each system may be found in section 5.1. The r|R e /ea 0 | 2 value 
was deduced from the absolute measurements (column 6 ) by use of the 
synthetic spectrum calculations described in section 3.5. These calcula- 
tions included all of the vibrational bands for which Franck-Condon 
factors are listed in appendix C. 
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TABLE E-l. TABULATION OF SHOCK-TUBE DATA FOR THE C 2 SWAN SYSTEM 
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TABLE E-2. TABULATION OF SHOCK-TUBE DATA FOR THE C 2 


DESLANCRES-d ' AZAMBUJA SYSTEM 
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L 


TABLE E-3. TABULATION OF SHOCK-TUBE DATA FOR THE C 2 PHILLIPS 
SYSTEM. Av = +2 SEQUENCE. 
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TABLE E-5. TABULATION OF SHOCK-TUBE DATA FOR THE C 2 BALLIK-RAMSAY SYSTEM 


V S » 

km/sec 

Q 

parts/cm 3 

W/cm z »u*sr 

E I Re/ ea o 1 ^ * 
A.U. 

m 

5520 

6.80*6 

59 

0.54 

i 

5550 

6.55*6 

77 

.65 

4.56 

5600 

6.00*6 

55 

.50 

4.70 

6090 

2.60*6 

37 

.71 

4.48 

5410 

8.05*6 

72 

.52 

4.45 

5400 

9.15*6 

91 

.90 

4.45 

5400 

9.15*6 

90 

.90 

4.35 

5205 

1.17*7 

113 

.90 

4.42 

5340 

9.90*6 

108 

.64 

4.32 

5160 

1.22*7 

133 

,65 

4.34 

5200 

1.18* 7 

188 

.94 

4.28 

5200 

2.00*7 

160 

.49 

4.35 

5315 

1.81*7 

260 

.83 

4.35 

5315 

1.81*7 

235 

.75 

4.40 

5405 

1.66* 7 

291 

1.00 

4.35 

5315 

1.81* 7 

258 

.82 

4.27 

5190 

2.01*7 

155 

.47 

4.36 

5330 

1.78*7 

235 

.76 

4.26 

5175 

2.03*7 

156 

.47 

4.38 

5365 

1.71*7 

196 

. 66 

4.34 

5300 

1.84*7 

223 

.71 

4.31 

5250 

1.91*7 

224 

.69 

4.37 

5350 

1.75* 7 

232 

.77 

4.49 

5495 

8.10*6 

65 

.46 

4.52 

5565 

7.20*6 

89 

.69 

4.51 

5540 

7.40*6 

93 

.72 

4.51 

5540 

7.40* 6 

71 

.55 

4.54 

5615 

6.60* 6 

81 

.69 
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TABLE E-6. TABULATION OF SHOCK-TUBE DATA FOR THE C 2 FOX-HERZBERG SYSTEM 


Test 

Pi. 

torr 

v.. 

km/ sec 

1 I 2 ’ 

K 

n C2 , 

parts/cnr 

I X‘ 

W/cm 2 *y *sr 

E|R e /ea 0 | 2 , 

A.U. 

1 

9. 

IS 

4.67 

5960 

3.25 16 

29.8 

0.38 

2 

9. 

15 

4.48 

5410 

8. 05 16 

43.8 

.46 

3 

10 

.0 

4.45 

5400 

9.15 16 

46.9 

.50 

4 



4.38 

5265 

1.1017 

48.7 

.36 

5 



4.42 

5340 

9.90 16 

53.5 

.38 

6 



4.45 

5400 

9. IS 16 

60.6 

.42 

7 



4.32 

5160 

1.22 17 

40.9 

.33 

8 



4.34 

5200 

1.18 17 

38.5 

.30 

9 



4.37 

5250 

1 . 1 2 1 7 

43.7 

.33 

10 



4.49 

5495 

8.10 16 

58.4 

.40 

11 



4.52 

5S65 

7.20*6 

62.0 

.42 

12 



4.51 

5540 

7.40 16 

68.2 

.46 

13 



4.51 

5540 

7.40 16 

67.2 

.45 

14 



4.54 

5615 

6.60 16 

66.9 

.45 

15 



4.52 

5565 

7.20 16 

50.5 

.34 

16 



4.54 

5615 

6.60 16 

46.0 

.31 

17 



4.44 

5390 

1 . 04 1 7 

67.8 

.49 

18 



4.47 

5445 

8.60 16 

55.6 

.39 


TABLE E-7. TABULATION OF SHOCK-TUBE DATA FOR THE Cz MULLIKEN AND 

FREYMARK SYSTEMS 


Test 

Pi, 

torr 

km/ sec 

Tz, 

°K 

NC2 / 3 
parts/ciy 

h> 

W/cm 2 »ysr 

E | R e /ea 0 | 2 , 
A.U. 


h- — 

Mul liken 

1 

10 

.0 

4.44 

5380 

1 . 04 1 7 

598 

0.17 

2 



4.49 

5495 

8.10 16 

407 

.12 

3 

1 


4.51 

5540 

7.40 16 

409 

.11 

4 

15 

.0 

4.39 

5380 

1.68 17 

705 

.12 

1 . Frevmark 
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A 

A *ir u 


uJl 


V 



b3 ^ 

C *1T 

g 

C ' \ 
c 

D 

? 

D h* 
u 

d 


E *E* 


E t U 


NOTATIONS 

constant appearing in equation (3.12), * h/8ir 2 cl^ 
lower electronic state of the Deslandres-d'Azambuja, 
Messerle-Krauss, and Freymark systems; upper electronic 
state of the Phillips system 

Einstein A coefficient for an electronic transition, 
sec -1 ‘particle -1 

Einstein A coefficient for the (v' , v") transition, 
sec -1 *particle -1 

lower electronic state of the Fox-Herzberg, Swan, and 

Ballik-Ramsay systems 

rotational constant, era -1 

black -body emission, W/cm 2 »ysr 

upper electronic state of the Ballik-Ramsay system 

upper electronic state of the Deslandres-d'Azambuja system 

upper electronic state of the Messerle-Krauss system 

speed of light, 2.9979 x 10 1 ® cm*sec -1 

rotational constant, cm -1 

dyadic of coordinate transformation 

upper electronic state of the Mulliken system 

total statistical weight or multiplicity 

upper electronic state of the Swan system 

energy, ev 

upper electronic state of the Freymark system 
integrated intensity in a rotational line, W/cm 3 *sr 
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e 


ea„ 



F(J) 


f e*« 


f v'v" 

G(v) 

G(X) 

5 

H 

H 0 

HO) 

h 


■h 



I 


J 

3 


spectral emissive power due to spontaneous emission, 
W/cm 3 *y*sr 

electronic charge, 4.80298 x 10" 10 dyne* cm 2 (esu) 
product of electronic charge and radius of first Bohr 
orbit, 2.5416 x 10" 18 esu*cm, unity in atomic units (A.U.) 
upper electronic state of the Fox-Herzberg system 
Gibbs free energy, ev 
rotational term energy, cm -1 

wavelength dependent electronic absorption f number, 
dimensionless 

band oscillator strength (absorption), dimensionless 

vibrational term energy, cm" 1 

radiometer calibration constant, V(W/cm 2 *sr) -1 

"gerade" or even state 

statistical weight of the level n 

Hamiltonian operator 

wavelength independent absolute strength parameter 
spectral irradiance incident on the radiometer, W/cm 2 *y 
Planck's constant, 6.623 x 10" 27 erg*sec, or enthalpy, 
cal/G 

Planck's constant divided by 2tt, 1.054 x io" 27 erg*sec 

heat of formation, kcal/mole 

moment of inertia, G*cm 2 

specific intensity, W/cm 2 *ysr 

total rotational quantum number 

total angular momentum 
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K 

It 

K 

k 

t 

L(X) 

A 

*s 

M 

«n 

m 

m e 

N 

ft 

P 

P 

Q 


* v" 

R 

ru£ 

IW? V , V ,.)I 2 

£|R«/«a 0 | 2 


rotational quantum number without spin 
total angular momentum apart from spin 
equilibrium constant 

Boltzmann constant, 1.38054 x 10“ 16 erg»°K -1 
orbital angular momentum 

spectral irradiance from the calibration source, W/cm 2 *p 

optical path length, cm 

radiating gas boundary, cm 

electric dipole moment, esu*cm 

mass of the nth nuclei, G 

particle mass, G 

mass of the electron, 9.1091 x 10~ 28 G 
particle concentration, particles/cm 3 
nuclear rotation vector 

pressure, atm, or rotational branch with AJ = -1 
parity, + or - 

partition function, dimensionless, or rotational branch 

with AJ = 0 

Franck-Condon factor 

rotational branch with AJ = +1 

transition moment, esu^cm 

electronic transition moment squared, esu*cm, 

3 *lMVv" )|2/du 

sum of the squares of the electronic transition moment, 
A.U., often referred to as the electronic transition moment 
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r 

r 


e 


r v'v" 

S 

-+• 

S 

Vv" 

SJ'A' 

J"A" 

qUV* J' 
*V"J" 

T 

T 

*e 

T(X) 

U 

U(r) 

u 


ee 


ne 
v nn 
V(r) 

V(X) 

v 

W g 

W, 


W,; 


internuclear separation, cm 

equilibrium internuclear distance (corresponds to the 

minimum of the potential curve), cm 
0 

r centroid, A 
spin quantum number 
spin angular momentum 
band strength, A.U. 

Honl-London factor normalized to (2S' + 1 ) (2J * + 1) 

total line strength, A.U. 

temperature, °K 

electronic term energy, cm -1 

normalized spectral response of a radiometer 

shock velocity, km/sec 

molecular potential function, ev 

"ungerade" or odd state 

volume, cm 3 

electron-electron Couiombic interaction, ev 

electron-nucleus Couiombic interaction, ev 

nucleus -nucleus Couiombic interaction, ev 

sum of the potential energy terms, * V * V + V , cv 
r 67 ee ne nn’ 

radiometer output, volts 
vibrational quantum number 

O 

Gaussian line width at half-height, A 

O 

Lorentzian line width at half-height, A 

O 

Voigt line width at half-height, A 
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Y 

Z 2 


°J 



0 

A 

t 

A 

*o 

V 

V 

e 

p 


lower electronic state of the Mulliken and Phillips systems 

dimensionless parameter, = A/B y 

dimensionless parameter defined by equation (3.81) 

rotational constant, cm" 1 

stoichiometric coefficients 

absorption coefficient, cm" 1 

absorption coefficient corrected for stimulated emission, 
cm" 1 

splitting constant, cm" 1 

Kronecker delta function, 5 = 1 for A = 0 and * 0 for 

O, A 

A/O 

splitting constant, cm" 1 
angle of oblique shock, deg 

quantum number of the angular momentum of the electrons 
along the intemuclear axis 

resultant orbital angular momentum of electrons along the 
intemuclear axis 

O . 

wavelength. A, or splitting constant (equation (3.43)), cm 1 

O 

band origin, A 

frequency of emitted radiation, sec" 1 
wavenumber, cm" 1 
wavenumber of bandhead, cm" 1 
wavenumber of band origin, cm" 1 
electron coordinates 
density, G/cra 3 


fa 
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E quantum number of the spin angular momentum along the inter 

nuclear axis 

E component of spin angular momentum along the internuclear 

axis 

o symmetry factor, =2.0 for homonuclear molecules, 1.0 for 

all others 

t radiative lifetime, sec -1 

<p lambda doubling function, cm -1 

y wave function 

0 quantum number of the resultant momentum along the 

internuclear axis 

total angular momentum along the internuclear axis, 

-+■ -*• 

= A ♦ E 



vibrational 

constant, 

cm -1 

0) X 

e e 

vibrational 

constant. 

cm" 1 

w y 
e 7 e 

vibrational 

constant, 

cm" 1 

w z 
e e 

vibrational 

constant , 

cm" 1 


Subscripts 

c calibration 

e electronic state or equilibrium position 

l lower state 

n nuclear state 

T test 

u upper state 


v 


vibrational level 
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f 

i 

f 

i 

I 

central point 

conditions ahead of shock wave 
conditions behind incident shock wave 
conditions in the driver tube 

Superscripts 

' upper state 

" lower state 

u upper state 

l lower state 


1 

2 

4 


f, 

X* 
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